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ABSTRACT 


The  desi^Vii  fatigue  resistant  structures  for  high  speed  aircraft  and  aerospace 
vehicles  depends  largely  on  the  prediction  of  realistic  acoustic,  fluctuating  aerody¬ 
namic,  and  engine  vibration  environments  and  on  the  estimation  of  the  attendant 
vibration  levels  of  structuiuil  components  and  attached  eejuipment.  The  practical 
engineering  limitations  on  the  mathematical  and  numerical  analyses  required  to 
treat  such  stiuctures  rigorously  by  classical  dynamics  necessitate  studies  of  alter¬ 
nate,  approximate  methods.  In  this  report,  a  definitive  statement  is  presented  of 
the  empirical  approacti  tni  detoi  iaining  ccrrclafieos  between  the  evr>;tntion  environ¬ 
ment  and  tlie  vibration  re.sponse  of  typical  fliglit  vehicle  structures  by  means  of 
statistical  analyses  of  measured  vibration  data.  The  various  aspects  of  the  vibra¬ 
tion  prediction  problem  and  tlie  general  philosophy  motivating  research  in  tlie  area 
of  empirical  correlation  arc  discussed.  Specific  treatment  is  given  to  the  effects 
of  bandwidth,  modal  density,  and  surface  pressure  spacecorni^lation  on  the  cross- 
correlation  of  energy  transmitted  along  various  structural  transmission  paths  in 
complex  linear  structures.  An  engineering  model  for  empirical  correlation  is 
developed  and  several  e.xainples  are  given  to  demonstrate  the  utility  of  the  equations. 
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SECTION  1 


INTRODUCTION 


At  the  present  time,  adequate  methods  do  not  exist  for  predicting  the  vibration  response 
levels  of  the  many  structural  and  equipment  components  of  aircraft  and  aerospace  flight 
vehicles.  Methods  generally  used  for  those  predictions  consist  of  vibration  analyses  of 
rather  elementary  structures,  such  as  uniform  beams,  plates  and  shells,  and  comparati’ 
analysis  in  which  measured  data  from  one  vehicle  are  used  to  estimate  response  levels 
of  another  vehicle  by  appropriate  mass  and  stiffness  scaling.  In  the  search  for  more 
powerful  prediction  techniques  and  more  realistic  vibration  response  estimates,  consid¬ 
eration  should  be  given  to  the  method  of  empirical  correlation,  in  which  statistical 
analysis  is  used  to  correlate  the  measured  responses  with  excitation  parameters,  overal 
structural  and  performance  parameters,  and  detailed  local  structural  parameters. 

The  purpose  of  this  report  is  to  present  the  general  philosophy  which  has  motivated 
research  in  the  area  of  empirical  correlations,  to  discuss  the  various  aspects  of  the 
vibration  prediction  problem,  to  develop  the  general  equa'ions  for  predicting  vibration 
responses  of  complex,  linear  structures,  and  to  set  forth  empirical  correlation  methods 
and  methods  for  developing  vibration  prediction  tools  from  these  correlations. 

Section  II  contams  a  general  discussion  of  the  engineering  problem  of  flight  vehicle 
structural  vibrations  and  the  relationship  of  this  problem  to  vehicle  reliabilitj'.  The 
environmental  prediction  problem,  the  vibration  respoiise  prediction  problem,  and  the 
importance  of  realistic  vibration  predictions  and  the  establishment  of  vibration  controls 
during  initial  design  arc  also  discussed  in  Section  II. 

In  Section  III,  the  general  equations  for  computing  vibration  :  espon.3e  in  any  complex, 
linear  structure  are  developed  for  point  force  excitation  and  distributed  pressure 
excitation,  having  discrete  and  continuous  frequency  spectra.  These  equations  are 
developed  in  terms  of  a  forcc-to-acccTcration  transfer  function  which  is  directly  relatabl 
to  transfer  impedances. 


Manuscript  released  by  the  .'luthors  July  19G4  for  publication  as  an  AFFDL 
Technical  Report. 
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Equations  for  moan- square.'  response  at  a  point,  autocorrelation  of  response,  cross¬ 
correlations  of  responses  at  two  points,  response  spectral  densities  and  cross- spectral 
densities  are  presented.  Si'veral  examples  arc  given  to  clarify  the  meanings  of  the 
transfer  functions  and  to  demonstrate  the  utility  of  the  equations. 

Some  essential  Ijaekground  material  on  the  development  of  empirical  correlation  tech¬ 
niques  is  presented  in  Section  IV,  Several  methods  for  obtaining  correlations  between 
excitation  and  response  using  statistical  techniques  are  developed  and  methods  for 
using  these  computed  correlations  as  vibration  prediction  tools  are  discussed. 


SECTION  II 


GENERAL  DISCUSSION 

Amons  the  many  technological  pi'o!)lcms  which  confront  the  aircraft  and  aerospace 
industries,  the  problems  associated  with  structural  and  equipment  vibrations  are  of 
paramount  importance.  Fatigue  of  structural  components,  and  malfunction  and 
breakage'  of  ''(luipment  in  flight  vehicles  often  accompany  high  vibration  levels;  such 
phenomena  can  lead  to  losses  of  mission  capability  and  even  to  catastrophic  failures. 
Past  and  present  generation  vehicles  have  experienced  severe  and  damaging  vibration 
response  to  engine  exhaust  noise,  engine  vibrations,  boundary  layer  turbulence,  and 
oscillating  shocks,  hi  several  instances,  time  consuming  engineering  studies  and 
costly  test  programs  were  required  to  reduce  the  levels  of  response  to  within  tolerable 
limits.  WUh  the  ever  increasing  sizes  of  vehicles  and  powerplants,  and  with  increasing 
flight  speeds,  the  vibration  excitation  environment  imposed  on  the  structure  may  becoim 
more  severe.  As  a  result,  the  vibi’ation  response  problem  for  futi  re  generation  flight 
vehicles  is  expected  to  become  increasingly  important  and  will  require  considerably 
greater  attention  than  has  been  devoted  to  this  problem  heretofore. 

VIBRATION  CQNTRO  L  riTROUGH  INITL^L  DESIGN 

The  operations  of  flight  vehicles  within  the  atmosphere  and  under  any  thrust  con¬ 
ditions  are  always  accompanied  by  sources  of  vibration  excitation.  Because  of  the 
flexibility  and  general  dynamic  nature  of  the  vehicle  structure,  a  portion  of  the 
oscillatory  energy  si!})plicd  by  the  excitation  sources  is  accepted  by  the  structure  and 
tramsmitted  thi  (jughout  the  vehicle.  Thus,  under  such  operational  conditions,  every 
structural  and  equiimient  component  is  placed  in  a  state  of  vibration  and  must  function 
properly  in  that  state.  Clearly,  vibrations  cannot  be  eliminated  completely;  hence, 
any  positive  apprtxu’h  towards  solving  the  vibration  problem  must  be  one  of  vibration 
control  instituted  for  the  purpose  of  minimizing  failure  and  malfunction  potentials 
and  thercljy  maximizing  the  overall  vehicle  reliability. 


Controlling  vibration  response  levels  by  controlling  the  environmentol  excitation 
levels  is  impractical  since  the  environments  aie  dependent  upon  such  major  design 
and  performance  parameters  as  engine  thrust,  engine  location,  vehicle  surface  con¬ 
tours,  trajectories  and  flight  speeds.  Rather,  vibration  controls  lie  in  the  design  of 
fatigue  resistant  structure,  the  design  of  equipment  to  withstond  operational  vibration 
levels  without  damage  and  malfunction,  the  positioning  and  mounting  of  equipment 
to  minimize  response  to  structurally  transmitted  vibratory  energy,  optimum  design 
of  structural  components,  the  reduction  of  stress  concentration  factors,  the 
development  of  maximum  damping  through  joint  design,  and  the  use  of  damping 
compounds.  As  such,  many  built-in  vibration  control  measures  can  be  incorporated 
doting  ll'.c  initio!  planning  and  dosig''  phases.  Other  detailed  vibration  control 
measures  can  be  used  during  the  advanced  design  and  qualification  testing  phases; 
while  still  oilier  limited  controls  can  be  supplied  after  vehicle  construction  from  a 
knowledge  of  vibration  response  levels  measured  during  operation. 

Although  the  need  for  vibration  control  and  the  extent  and  degree  of  the  control 
cannot  bo  ptedel.ciminc'd  fer  all  key  structural  ond  equipment  components,  the  most 
effective  and  the  most  economic  controls  are  those  Incoiporoted  in  the  initial  de¬ 
sign  of  both  the  structure  and  the  equipment.  Altering  the  structure  after  fobrication 
begins  is  difficult  c;nd  in  many  Instances  it  is  hard  to  justify  on  the  basis  of  the  value 
gained  when  weighed  ogoinst  the  time  and  expense  required  for  changing  drowings, 
toolinp,  and  comUnietion  pr^'ceduro.s.  The  usual  result  is  a  compromise  in  which  the 
full  effect  of  the  newly  Instituted  control  measure  cannot  be  realized.  This  is  gener¬ 
ally  less  true  of  cc]jlpment  packages  which  are  not  integral  parts  of  the  structure  and 
which  can  be  easily  removed;  however,  here  again,  valuoble  time  may  be  required 
for  cvaluntiun  of  the  viljration  problem,  redesign,  and  retesting  for  tjualification. 
Consideration  of  the  vibration  response  problem  mrrst  therefore  begin  during  prelimi¬ 
nary  design  . 

Ideally,  these  vonsiderutlons  should  begin  with  the  estimations  of  environments,  at 
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the  vehicle  surface,  caused  by  such  fluctuating  pressure  sources  as  engine  exhaust 
noise,  boundary  layer  turbulence,  separated  flow,  base  pressure  fluctuations,  jet 
wake  noise  and  oscillating  shock  waves.  The  various  characteristics  of  Interest  for 
these  environments  are  distributions  of  pressure  levels  over  the  vehicle  surface, 
frequency  distributions,  propagation  speeds,  impingement  angles  and  spatial  corre¬ 
lations.  Estimates  should  also  be  made  of  the  direct  vibratory  input  to  the  structure 
due  to  engine  burning,  vibrations  of  propellant  piping  system',  and  other  on-board 
mechanical  sources  of  excitation .  From  the  knowledge  of  the  estimated  environmen¬ 
tal  I  evelsand  other  environmental  charocteristics,  vibration  response  levels  would 
be  predicted  for  vibration  sensitive  structural  components  In  regions  of  high  level 
excitation  and  for  all  key  structural  and  equipment  components  such  as  skin  panels, 
skin  stiffeners,  brackets,  piping  systems,  equipment  mounting  structure  and  equip- 
rnerit  comporents.  Specific  requiiements  for  vibration  control  could  then  be  defined 
from  the  estimated  response  levels,  structural  fatigue  criteria,  and  equipment  mal¬ 
function  criteria.  Finally,  detailed  design  techniques  could  be  developed  In  order 
to  obtain  sufficient  vibration  control  to  reduce  the  severity  of  potentially  damaging 
vibrations;  and,  realistic  qualification  test  specificotions  could  be  established  for 
critical  structural  and  equipment  components.  Such  an  effort  should  result  in  the 
maximum,  practical  built-in  vibration  control.  Continuing  this  effort  through 
the  various  refinement  phases  of  advanced  design,  qual 'f  ication  testing  and  flight 
testing  should  result  in  the  maximum  vehicle  reliability  within  the  limits  that  vi¬ 
bration  induced  failures  influence  overall  vehicle  reliability. 

The  three  principal  tools  required  for  attaining  vlb'^ation  control  through  initial  de¬ 
sign  are  methods  for  predicting  the  environmental  levels  ond  other  environmental 
parameters,  methods  for  predicting  vibrotion  response  levels  and  transmission  char¬ 
acteristics  of  vehicle  structures,  and  design  techniques  for  minimizing  structural 
response  levels  and  safeguarding  equipment  against  damage  and  malfunction.  During 
the  past  ten  years,  in  which  significant  odvances  have  been  made  in  the  development 
of  large  turbojet  and  rocket  engines  ond  in  attaining  higher  flight  speeds,  a  consider- 
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able  engineering  effort  has  been  devoted  to  the  measurement  and  analysis  of  exhaust 
noise,  boundary  layer  turbulence,  and  other  fluctuating  pressure  phenomena.  From 
tills  effort  have  come  empirical  prediction  methods  which  allow  for  reasonably  accurate 
estimates  of  tlie  environmental  levels,  distribution  over  the  structure,  frequency  dis¬ 
tributions,  propagation  characteristics,  and  in  some  cases,  spatial  correlation  properties. 

In  contrast,  practical  engineering  methods  for  predicting  vibration  response  levels 
in  typical  flight  veliicle  structure  have  not  been  developed  primarily  because  of  the 
complexities  of  the  stioictural  acceptance  and  transmission  properties.  The  few 
limited  attempts  to  develop  empirical  vibration  prediction  methods,  based  on  measured 
data,  show  several  orders  of  magnitude  of  scatter  in  the  data  even  when  normalized 
against  typical  weight,  stiffness,  and  excitation  parameters,  A  typical  scatter  in  the 
vibration  data  may  be  20  decibels  whereas  a  tj’pical  scatter  in  appropriately  norma¬ 
lized  environmental  data  may  bo  2  to  .3  decibels.  The  sensitivity  of  the  vibration 
levels  to  local  structural  parameters,  which  can  have  wide  variations  throughout  a 
structure,  and  to  the  highly  resonant  character  of  a  structure,  is  considered  to  be  a 
major  factor  in  causing  such  a  broad  scatter.  Generally,  response  predictions  are 
made  by  idealiziag  complex  structural  components  to  mathematically  tractable 
systems  and  analy/ing  such  systems  by  cither  lumped  parameter  methods  or  by  general¬ 
ized  harmonic  analysis  methods.  Except  for  the  use  of  damping  compounds  and  the 
design  of  joints  for  high  damping,  most  design  techniques  for  vibration  control  are 
dependent  upon  a  knowledge  of  the  dynamic  characteristics  of  the  structure.  A 
major  effort  must  therefore  be  made  to  dcwclop  vibration  response  prediction  methods 
which  will  account  for  necessary  excitation  parameters  and  major  and  local  structural 
parameters. 


THE  MAJOR  FREQUENCY  REGIMES  OF  STRUCTURAL  RESPONSE 

The  Low  Frccniency  Ui'^ime  -  Overall  Fundamental  Vehicle  Modes 

The  fundamental  principh'S  of  elasticity  and  mechanics  provide  alt  of  the  analytical 
relations  that  are  luajuireil  to  mathematically  describe  the  dynamic  response  of  any 
structure,  sinipdi'  or  eomi>lex,  by  means  of  tlifferential  equations  of  motion.  However, 
since  the  dynamic  al  complexity  of  llu'  structure  has  a  definite  bearing  on  the  complexity 
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of  tho  equations  of  moMon  and  their  solutions,  a  direct  mathematical  approach  to  the 
vibration  prediction  problem  is  limited  to  those  systems  which  display  more  or  less 
elementary  rest  nance  characteristics.  In  the  sense  used  here,  elementary  resonance 
characteristics  refer  to  resonant  deflection  shapes,  or  mode  shapes,  which  can  be  ex¬ 
pressed  mathematically  in  closed  form  or  which  can  be  obtained  with  reasonable  accu¬ 
racy  by  numerical  techniques,  such  as  finite  difference  methods,  or  by  lumped  param¬ 
eter  teclmiques.  Examples  of  such  systems  for  which  the  mode  shapes  can  generally 
be  expressed  in  closed  form  are  uniform  rods  with  longitudinal  resonances,  uniform 
beams  and  plate  i  'vith  bending  resonances,  and  uniform  shells  with  uncoupled  exten- 
sional  and  bending  resonances. 

Other  systems  wliich  are  complex  in  geometry,  mass  distribution  or  stiffness  distribu¬ 
tion  can  have  elementary  resonance  characteristics  for  the  first  few  resonant  modes. 
Although  the  mode  shapes  will  seldom  admit  to  simple  mathematical  represcnUition, 
numerical  techniques  and  lumped  parameter  methods  can  be  used  to  easily  define  the 
motie  shapes  and  compute  llic  resonant  frequencies.  E.xamples  of  such  systems  are 
overall  vehicles  vibrating  in  the  fundamental  bending  mode  or  fundamental  longitudinal 
mode,  and  stiffened  panel  configurations  vibrating  in  overall  configuration  modes  whose 
resonant  frcquc'ncies  are  well  below  those  of  the  individual  subpanels.  Ihe  significant 
feature  to  be  observed  in  the  latter  two  e.xamples  is  that  the  elastic  wave  lengths  at 
resonance  are  long  relative  to  the  dimensions  and,  hence,  to  the  resonant  elastic 
wave  lengths  of  the  individual  structural  components,  Conscfjuenlly,  the  complex 
local  variation.s  in  mass  and  stiffnc'ss  are  averaged  out  so  that  much  of  this  fine  detail 
can  be  neglected,  'riuis,  lor  c.xample,  a  panel  having  a  large  number  of  closely 
.spaced  orthogonal  stiffeners  can  be  treated  as  an  equivalent  uniform  panel  in  the  funda¬ 
mental  jnode.  Higher  modes  of  the  panel  may  be  similarly  approximated  with  the  restri 
that  the  overall  stiffened  [)an(d  wave  lengths  are  much  longer  than  the  distances  between 
stiffene  rs. 

T ild.! rj'i’iElfiniil Hang(‘  -  Fundamental  Resonances  of  Stnictur.al  Elements 

When  l)uilt'up  structures  are  e.xcited  at  frequencies  above  those  of  the  overall  resonant 
modt'S,  the  (da.^dic  biuvling  aivl  torsional  wave  lens-dhs  are  of  the  same  order  as  the 


dimcMisions  cif  individiml  sti'uctural  componcMits.  At  these  wave  lengths  the  local 
variations  in  stiffness  and  mass  distribution  significantly  affect  the  local  character 
of  the  vibration  response.  As  a  result,  the  higher  order  resonances  of  the  vehicle 
tend  to  breal-  up  into  resonances  of  individual  struciural  components,  or  small  groups 
of  such  components,  so  that  the  overall  modal  pattern  becomes  exceedingly  complex. 

It  is  impractical  at  these  frequencies  to  consider  the  local  resonant  modes  as  being 
an  integral  part  of  some  liigher  order  resonant  mmles  of  the  overall  vehicle.  In  fact, 
two  structural  components  which  arc  wididy  separated  may  be  essemtially  uncoupled 
so  that  their  responses  are  nearly  independent.  Present  methods  for  predicting 
vibration  response  in  this  frecpieney  range  generally  consist  of  the  analysis  of  the 
response  of  single  eompononts  to  the  local  tmyironment  acting  on  the  eomjionents. 
Tj’pically,  tlu'  structural  components  aii'  treated  as  uniform  beams,  plates  and 
shells  for  which  matluonatica!  solutions  of  the  i-esponse  to  fluctuating  surface  pressures 
or  to  other  fluctuating  excitations  are  known. 

The  accuracy  of  such  nn'tliod.s  depends  upon  the  aeeuraev  with  which  tiie  exeitatioii 
can  be  mathematically  slniulateil.  the  degree  of  accuracy  in  approximating  l!ie  tnu' 
natural  boundary  conditions  of  tlic  component,  and  the  degree  of  dynamical  inde¬ 
pendence  of  the  component  with  re.speet  to  neighboring  components.  .Structural 
components  such  as  skin  panels  excited  by  fluctuating  surface  prcjvsures  arc  sensitive 
to  the  ratios  of  elastic  wavt  lengHi  ti'  pressure  wave  lengtlis.  pressure  projnagation 
speeds,  angles  of  impbigemcnt  and  spatial  correlation  lengtlis.  These  e.xeitatiun 
parameters  must,  tl'.crefore,  be  Icnowii  in  order  to  havi*  reason;i])ly  accurate  response 
predictions. 

The  re.sonant  freqiieneies,  detlection  anqilitudes,  and  strain  levels  of  panels  are  also 
dependent  upon  tl’c  lioundary  condition.s  at  the  edges  of  the  panel.  The  bending 
and  torsional  stiffnesses  of  the  edge  memtiers  which  support  the  jianel  determine  tlie 
degree  of  translational  and  rotational  fixity  of  tlie  panel,  and  the  edge  deflections 
associated  with  the.se  flexibilities  have  a  marked  effect  ujion  strain  amplitudes  at  llu' 
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edge.  Tlic  dynamical  properties  of  these  edge  members  also  control  the  flow  of 
vibratory  energy  from  the  panel  to  neighboring  structural  components  and  from  the 
neighboring  compoiumts  to  the  panel.  Wlic  i  the  fundamental  bending  and  torsional 
frequencies  of  these  edge  stiffeners  are  much  higher  than  the  fundamental  panel 
frequency,  response  of  the  first  few  modes  of  the  panel  can  be  predicted  fairly 
accurately  by  treating  the  panel  as  an  individual  component  which  is  essentially 
uncoupled  from  the  surrounding  structure.  However,  if  the  flexibilities  of  the 
edge  members  are  significant,  a  group  of  several  panels  may  have  to  be  analyzed 
simultaneously.  This  is  particularly  true  in  the  case  of  regular  structure  in  which 
several  connecteti  [lancls  have  nearly  the  same  resonant  fretiuencies  --  the 
individual  structural  components  do  not  respond  independently  but  are  highly  coupled 
with  adjace  nt  components. 

Those  consideration.s  are  not  confined  to  skin  panels  but  are  valid  for  any  single 
structural  component  or  group  of  components  within  tlie  vehicle'.  In  general  terms, 
if  the  supjiort  structure  and  structural  components  adjacent  to  a  certain  structural 
component  e.xhii 't  a  high  impedance  relative  to  the  impedance  of  the  component,  then 
the  response  of  the  component  will  be  essentially  independent  of  the  n.cigliboring 
structure,  Wlicn  the.se  impeclanees  are  of  the  same  order  of  magnitude,  dynamic 
coupling  can  be  expected. 

Although  many  vibration  analyses  in  the  medium  frcciuency  range  are  bused  on  such 
idealizations  as  uniform  Iieams  and  plates  having  pinned  anri  clamped  edges,  serious 
efforts  are  being  made  to  develop  more  general  methods  which  account  for  the  effects 
of  local  jiaramctcrs  and  tiius  allow  for  oetter  appro.ximations  of  the  natural  Imundary 
conditions  and  tlie  dynamic  coupling  of  adjacent  structural  components.  It  is  difficult 
to  predict  how  far  tlic  .‘^tatc  of  the  art  of  viliration  ana’ysis  of  complex  .structures  cuiibe 
exteniiod  by  if.eans  of  such  ga  niM-alizat ions.  The  increasing  storage  capacity  and 
c  omputing  speeds  of  digital  compu.ers  allows  the  dynamicist  to  inciude  an  ever  larg  r 
numlier  of  strut  tui'al  mass,  stiffness  and  geometry  parameU'rs  and,  in  time,  quite 
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complex  systems  may  be  analy/.ahic.  Two  dimensional  stiffened  panel  and  shell 
configurations  will  certairly  be  amendable  to  numerical  computation  in  reasonable 
computing  times  with  future  generation  computers.  Bulkhead  and  some  internal 
structure  will  probably  also  be  included  at  allowable  computing  costs. 

Unfortunately,  the  number  of  parameters  and  the  general  size  of  the  computing 
task  increases  very  rapidly  with  the  number  of  individual  structural  components 
being  considered  and  with  the  order  of  the  response  mode  being  analyzed.  There 
will  be  at  least  five  parameters  reijuired  to  determine  the  dynamics  of  each 
structural  element  included  in  the  analysis,  and  several  resonant  frequencies  will 
be  required  for  each  such  element.  As  the  order  of  the  resonance  increases, 
the  amount  of  computing  required  to  locate  system  resonant  frequencies  increases 
because  of  the  nature  of  tlie  numerical  methods  used  for  finding  these  resonances. 

In  addition,  the  modal  density  incnuises  with  frequency  ai  that  the  resonant 
frequencies  become  closely  spaced  and  mcxlal  deflection  patterns  become  increasingly 
hard  to  display.  Complexity  in  the  modal  patterns  of  the  skin  panels  implies  difficult 
and  lengthy  computations  of  the  griienilized  or  'effective  force  caused  by  random 
fluctuating  surface  pressures.  Thus,  although  the  computer  techniques  are  expected 
to  provide  a  bright  future  for  the  analyst,  it  is  clear  that  there  exists  a  definite 
upper  limit  on  the  degree  of  complexity  that  can  be  ii.cluded  in  vibration  analyses  for 
the  medium  frequency  range. 

The  High  Frcciuency  Regime  -  Modal  Response  in  Smaller  Structural  Elements 

At  frequencies  higher  than  the  first  few  resonant  modes  of  individual  structural 
components,  the  modal  density  at  all  points  in  a  complex  vehicle  structure  is  so 
higji  as  to  warrant  any  exact  analysis  as  impractical.  In  this  frequency  range,  com¬ 
plex  coupled  resonances  will  occur  in  skin  panels,  panel  stiffeners,  stiffener  flanges 
and  webs,  bracketry,  equipment,  bulkheads,  piping,  etc.  As  such,  the  structural 
response  is  similar  to  a  reverberant  acoustic  field  in  that  vibratory  energy  can  be 
easily  transmitted  in  all  directions  throughout  the  structure.  In  part,  the  difference 
between  the  transmission  properties  in  tlie  high  frequency  region  and  those  in  the 
mf’dium  frequency  range  i  ^  due  to  the  fact  that  the  natural  stiffness  barriers  which 
tend  to  confine  the  vibratory  energy  to  a  single  component,  such  as  a  skin  panel,  or 
a  group  of  structural  components,  no  longer  are  effective.  Those  members  which 


provided  the  stiffness  barriers  are  themselves  in  resonance  and  can  thus  easily 
transmit  energy  at  these  frequencies.  Predicting  vibratory  response  levels  in 
such  a  frequency  range  might  more  readily  be  accomplished  by  considering  the 
structure  as  a  transmission  line  and  treating  tlie  vibration  response  as  a  random 
process  to  be  described  in  terms  of  statistical  parameters. 

CONSPECTUS  OF  THE  STATISTICAL  APPROACH  TO  VIBRATION  PREDICTION 


Analytical  techniques  involving  the  solutions  of  equations  of  motion  of  a  structural 
system  provide  the  necessary  vibration  prediction  tools  for  the  low  frequency 
regime  of  the  structure  and  for  certain  mathematically  tractable  structural  con¬ 
figurations  in  the  medium  frequency  regime.  For  highly  coupled,  complex 
configurations  in  the  medium  frequency  regime  and  for  high  frequency  excitation 
of  the  structure,  other  less  exact  methods  must  be  developed  for  vibration  pre¬ 
diction.  The  use  of  statistical  techniques  and  the  development  of  statistical 
relations  between  vibration  response  levels  and  the  significant  parameters  which 
describe  the  excitation  and  the  structure  appears  to  be  the  only  reasonable  approach 
to  so  complex  a  problem.  Such  relationships  should  be  developable  in  terms  of 
correlation  coefficients  and  other  statistical  averages  from  a  knowledge  of  measured 
vibration  data,  measured  or  estimated  environmental  data,  and  from  a  knowledge 
of  the  structural  configuration  on  which  the  responses  were  obtained.  Once  these 
relationships  are  developed  and  the  necessary  numerieal  evaluations  are  made  for 
representative  flight  vehicle  structure,  a  practical  engineering  method  should  result 
for  vibration  response  prediction  during  initial  design. 

In  order  to  develop  the  appropriate  statistical  relations  it  is  necessary  first  to  identify 
all  of  the  significant  excitation  and  structural  parameters  that  might  significantly 
influence  vibration  response.  These  parameters  might  include  engine  thrust,  flight 
speed,  altitude,  mass  of  vehicle,  vehicle  surface  area,  average  size  of  skin  panels, 
distance  from  engines,  distance  from  projections  and  aerodynamic  discontinuities, 
frequency,  pressure  wave  lengths,  length  of  vehicle,  skin  panel  thickness,  component 
mass,  component  stiffness,  distances  between  stiffeners,  etc.  The  actual  relations 
between  these  parameters  and  response  amplitudes  can  then  be  estimated  by 
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Sf'voral  (itlTerc'U  nielhj  ls,  once  a  reijrosenlative  sampling?  af  mcasure'l  vibration 
data  is  available.  Using  linear  regression  anjilysis  techniques,  linear  relationships 
can  be  conipulcJ,  on  a  statistical  basis,  between  response  amplitudi’  and  each  of  the 
above  excitation  and  structural  parameters  taken  singly  or  as  a  group.  Or,  definite 
mathematical  relations  could  ae  established  between  the  response  amplitude  and 
•.hese  pai'ameters,  and  a  statistical  amilysis  can  be  conducled  to  determine  how  well 
the  measured  data  fits  the  mathematical  model.  In  the  latter  case  the  closer  the 
mathematical  model  de.scribes  the  vibration  phenonu  na,  the  closer  will  be  the  fit. 

As  another  alternative,  very  general  but  realistic  mathematical  relations  could 
be  established  containing  unknown  constants  which  can  be  determined  so  as  to 
obtain  a  best  fit  between  the  relationships  and  ’he  major  body  r  f  the  measured 
data.  If  a  broad  and  representative  sampling  of  measured  vibration  data  could  be 
accumulated,  and  if  these  throe  methods  could  be  used  for  various  types  of  structures 
(such  as  fighter  aircraft,  bombers,  small  missiles,  large  rocket  boosters)  and  for 
various  flight  regimes,  statistica.1  relationships  could  be  developed  which  would  ue  of 
significant  value  as  vibration  p^'cliction  tools. 
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SECTION  III 


STRUCTURAL  VIBRATION  RESPONSE  THEORY 


In  developing  empirical  correlation  techniques  for  analyzing  the  vibration  responses 
of  complex  structures,  it  is  necessary  to  obtain  ex})ressions  for  the  response 
amplitudes  of  the  structure  in  terms  of  the  various  structural  and  excitation  param¬ 
eters  that  influence  vibration  response.  For  example,  it  is  necessary  to  have 
functional  relationships  between  the  amplitude  of  acceleration  response  at  a  point 
within  the  structure  and  the  frequency  of  excitation,  amplitudes  of  fluctuating 
acoustic  and  aerodynamic  pressures,  distributions  and  spatial  crosscorrelations  of 
these  pressures  on  the  surface  of  the  structure,  input  impedances  over  the  areas  of 
excitation,  and  transfer  impedances  between  the  excitation  and  the  point  within  the 
structure  where  response  occurs. 

Such  equations  can  bo  obtained  by  several  different  methods  and  can  be  w'ritten  in 
several  different  forms.  Partial  differential  equations  that  express  the  exact 
balance  of  elastic,  inertial,  and  applied  forces  can  bo  set  up  for  each  structural 
element  and  the  response  of  any  element  can  be  treated  ns  a  boundary  value  problen 
With  this  method  it  is  necessary  to  equate  deflections,  slopes,  bonding  moments  anc 
shears  of  adjacent  structural  elements  along  the  lines  of  attachment.  Alternatively, 
the  structural  response  can  be  expanded  in  terms  of  the  normal  modes  of  the  structi 
and  ordinary  differential  equations  can  bo  solved  for  the  amplitude  of  response  of 
each  mode.  In  this  method,  it  is  necessary  to  determine  the  individual  mode  sh.apcf 
the  generalized  mass,  generalized  damping,  and  generalized  force  for  each  mode. 
However,  in  the  analysis  of  complex  structures  having  complex  random  forcing  func 
it  is  necessary  that  the  final  equations  to  be  used  contain  only  quantities  that  are 
known,  that  can  be  estimated  with  reasonable  aceuraev,  that  can  be  measiu’ed  on 
a  practical  engineering  basis,  or  that  can  be  determined  statistically  from  measure' 
vibration  and  excitation  data.  Since  the  numerous  local  structura’  details  and  the 
exact  transmission  characteristics  of  a  complex  structure  are  generally  unknown,  i 
appeal's  that  the  equations  of  interest  must  be  e.xpressed  in  terms  of  input  and  trans 
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imped.'incc.s  or  other  similar  quantities  which  account  for  the  collective  effect  of 
mass,  stiffiK'ss,  and  dampinj^  distributions  between  points  in  the  structure.  For  the 
pi'oblem  at  hand,  equations  of  this  kind  are  more  appropriate  than  thf)se  involving 
the  e.xac't  balance  of  elastic,  inertial,  and  applied  forces  in  the  structure,  or  those 
which  require  a  knowlca4;e  of  the  very  complex  mode  shapes  of  the  vehicle  in  the 
medium  to  high  frequency  I'ange. 

In  this  section,  the  acceleration  response  of  an  arbitrary,  linear,  deterministic 
structure  is  an.aJyzed  for  various  tjTios  of  loading.  Expressions  are  derived  for  the 
autocorrelation,  spectral  density,  and  mean-square  acceleration  response  at  a 
point  on  *he  structure.  Ex})ressions  are  also  derived  for  the  crosseorrelations  of 
applied  forc'C  and  ''esi)onse  and  for  the  crosseorrelations  of  the  acceleration  rc.sponses 
between  point:  on  the  structure.  These  e.xiiressions  are  developed  for  both  multiple 
point  force  loadings  and  distributed  pressure  loadings.  In  eac-h  ease,  discrete  and 
continuous  frequency  spectra  arc  considercti. 

The  structure  is  a.‘<.sunu'd  to  be  linear  in  order  to  lake  a<lv:^n1age  of  the  principle  of 
superposition  of  rcspon.s(\s  at  a  point  due  to  several  vibratory  inputs  acting  over 
different  parts  of  the  structure.  In  the  ease  of  nonlinear  structures,  the  linear 
superposition  principle  is  invalid  and  special  methods  must  be  used  to  treat  each 
different  Ivpe  of  nonlinearity.  Nonlinearities  are  e.xpected  to  occur  in  certain 
limited  regions  of  the  structure;  however,  on  the  average  over  the  whole  structure, 
linearity  is  expected  to  dominate.  The  structure  is  deterministic  in  the  sense 
that  all  input  impedances,  transfer  impedances,  and  other  quantities  appearing  in 
the  response  e.xpi’essions  arc  assumed  to  be  known.  In  this  sense,  the  expressions 
to  be  developed  are  exact  for  linear  systems.  A  discussion  of  the  necessaiy 
simplifications  and  approximations  required  to  statistically  analyze  vibration 
re.spon.‘'0  of  complex  structures  is  reserved  for  Section  IV. 


H 


POINT  FORCE  EXCITATION  AT  DISCRETE  FRECUENCIE5 


Force-to-Acceleration  Transfer  Function 


The  net  vibration  ttJsponso  at  any  point  in  a  structure  is  the  cumulative  effect  of  all 
of  the  oscillatory  forces  acting  on  the  structcie.  Forces  applied  at  distinct  points 
induce  local  vibrations  at  the  points  of  appi icat ion;  and  these  local  vibrations 
are  transmitted  along  common  structural  paths  and  add  to  give  net  responses  in 
accordance  with  their  relative  phases.  Vvlmn  the  response  amplitudes  are  linearly 
dependent  upon  the  amplitudes  of  the  applied  forces,  the  response  amplitude  at  a 
point  Is  the  linear  summation  of  the  responses  at  that  point  which  result  from  ex¬ 
citations  at  each  pc’nt  on  the  stracture. 

Consider  for  example  the  structural  system  shown  its  Figure  1  which  has  forces  F^(t) 
applied  at  points  sT  ,  and  which  has  o  structure!  element  at  whose  net  acceleration 
response  is  V,Mr).  It  is  assumed  that  coch  of  the  point  forces,  ^^^0,  is  cr  steady- 
state  force  having  several  discrete  frequency  components,  frequencies 

f|^  ,  so  that 


F  (t)  F  ,  (tt  (I) 

n  /  j  nK 

h 

If  denotes  the  component  of  W^(t)  resulting  from  the  force  component  F 

then  for  a  linear  structure,  the  net  acceleration  re'r^onsc  of  the  element  ot  y  is 

r 

given  by  the  expression 


E  E 

k  n 

where  the  summation  Index  n  extends  over  all  applied  forces  and  where  the  summation 
index  k  extends  over  oil  discrete  frequencies,  f^^  . 
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The  relolionsliip  bctwcrn  tlie?  applied  foice  and  the  vibration  response  caused  by  this 
force  involves  the  or  cenfoiu.r  ond  transmis'ion  characteristics  of  the  structure  between 
the  points  of  loud  npplicntiofi  rinri  lesponsi'  .  Because  of  tlie  sensitivity  of  the  structure 
to  excitotlijn  f  1  equf'nc i es^  the  stiuctuie  its  ts  os  o  complex  filter  which  accepts  and 
transmits  rertain  frequencies  more  readily  than  others;  and  hence  the  relationships 
between  the  applied  forces  and  the  responses  ore  dependent  upon  trie  excitation 
frequency.  Once  these  r e lot Icjnshl ps  are  known  for  each  frequency  component  of 
each  of  tire  np()lled  lotre'  ,  the  net  acceleration  resfxsnse  at  a  qiven  point  can 
be  expressed  in  terrrrs  of  th<'  amplitudes  and  relative  phases  of  the  individual  forces 
and  in  terms  of  thf'  sturc  turnl  oc ceptacires  and  transmission  characteristics  from  each 
inpt.it  prslnt,  ,  to  the  tjiven  tesuonse  point  y  .  The  steody -state  excitation, 

ti  t 

F  ,  (t)  ond  the  steudy-state  response  w  ,  ft)  c:at>  be  related  by  inttodircinq  o  foice- 
nv  ntr 

to-occeleration  transfer  frrrrctlorr,  a  (f  )  which  accoijnts  for  tire  irrprjf  and  tronsfer 

nr  k 

imrsedorrc.r  properties  of  the  structure  at  fteqtiency  f,  .  Irr  terms  of  o  (f,  )  /  the 

'  ^  rr  r  l< 

excitatlrrn  trrnrl  respotr'e  cue  related  by  the  eq'jcrtlon 


w  (tl  a  ff,  '»  F  ,  ftl 
ntk  nr  v  nk 


(3) 


The  net  acceletatiorr,  'A'  ftt,  of  the  sfrurrotal  element  crt  y  is  ohtnined  by  stih- 

I  I 

stitutirrf;  Eqrroflon  13)  inter  EqriOtlon  (?)  qi  vin'.) 

vi  (t)  V'  Cl  ff,  I  F  ,  ft)  f4) 

I  ^  nr  k  nk 

k  n 


lire  fprontity  u  ff,  (  ,  which  Ins  the  units  of  inverse  mass,  is  related  tc)  the  stri.rctural 
rtr  k 

troo'-fer  i 'irpedonce ,  7  ff  )  ,  by  the  expression 

nr  k 


a  ft,  ) 
nr  k 


2n  I  t. 


/  0.) 

nr  k. 


(5) 
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where  ttie  tronsfer  impedance  ot  frequency  ic-  defined  by  the  equation 


F  Jf)  Z  (f.  )  w  Jt) 
nk  in  k  nrk 


(6) 


Like  the  transfer  impedance,  t!ie  quantity  a  (f^^)  is  complex,  and  can  be  written 


tlie  form 


a  (f.  )  la  (f,  )  I  0 
nr  k  I  nr  k 


i0  (f,  ) 
nr  k 


(7) 


wfiete  o  (f.  ^  Is  the  modulus  and  ft  (f,  )  is  the  phase  anqie  of  a  (f,  ) 
nr  k  nr  k  ni  k 


To  illustrate  the  functional  fotm  iif  u  (f  )  foi  simtrle  i/'-iems,  consldei  the  one- 

ni  k 

deqree -of-fieedom  spi  i  n;) -mns'.  system  shown  in  Figure  2.  Let  the  applied  force, 

F ft),  be  stcady-siote  and  harmonic  at  the  single  frequcricy  f  .  Then  os  shown  on 
page  1-7  of  Reft  rence  1,  the  steady-state  di'flection  irrsponse  V/(t)  of  the  moss  i: 
fjiven  by  the  Equation 


r(n/k 


whei  e 


and  I, 


(1/  Ztt)  ^ 


:/  m 


Thus  the  acceleration  response  V/ft)  t  nn  be  wiitten  in  tlie  form 


(3) 


in 


17 


Compoiing  Eajnfioii  (8)  witi>  Eq^Jtion  (4:  v.viv.-.  »8o-t  Uie  !>onsfei  function  off)  for 
fhf-  limpU'  spt  ing -mo'.',  '  .itn-n  i5 


o(n  -  H(f )  , 

witerc  the  sinqie  dogrf^t; -»:f-freoriorn  resfy^ns''  ♦s.incfion  H(fl  ii 


(9) 


The  modulus 


and  the  phase  angle  l^ff)  of  off)  ore 


(10) 


|a(f)|  (?T,f)'^  |H(f)|,h 


fliequontity  jhlffljin  Equation  (11)  is  the  dynoTiir  moqni  ticatlon  factor  for  the 
single  degree -of -freedom  system,  and  this  function  fios  o  maximum  value  of  O  at 
resonance.  As  shown  in  Reference  1,  poge  1  i'n  function  H(f)  is  the  Fourier 
transform  of  the  unit  Impulse  functiori  hft),  that  i' 

Hff)  k 

-a' 

where  tite  unit  impulse  function  for  the  single  .jopifca  -if  tiopd  system  Is 
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1  ‘^V 

h(t)  - -  e  sint  t  f  0 

m  <T 

-  0  t  40 


(13) 


The  quantity  f  Is  the  damping  ratio,  a  Is  the  damped  natural  frequency,  and  Uq  is 
the  undamped  natural  frequency,  where 


f -  1/2G  , 


Mg. 


v; 


1 


Thus  from  Equations  (9)  and  (12)  it  is  seen  that  the  transfer  ’’unction  a(f)  can  be 
written  in  terms  of  the  Fourier  transform  of  the  unit  impulse  function 


a(f)  -  (2i»f)^ 


7 

-00 


,  ,  .  -2nift  , 
h(t)  e  dt 


To  illustrate  the  functional  form  of  a  (f.  )  for  a  second  example,  consider  the 

nr  K 

uniform  bv.  m  shown  in  Figure  3  which  has  a  single  harmonic  force  F(t),  of  frequency 
f,acting  at  point  x  .  The  deflection  response  V/ft)  at  point  y  on  the  beam  con  be 
found  from  the  normal  mode  expansion 


W(t)  ^  q^  (t)  (y)  (14) 

I’ 

where  <p^(y)  is  the  >■ -th  mode  shape,  normalized  to  unity,  and  q^,  (t)  Is  the  deflection 
amplitude  of  the  r  -th  mode.  V/hen  the  modes  respond  independently,  the  generalized 
coordinate  q^  (t)  of  the  beam  satisfies  the  single  degree-of  freedom  equation 

w,.  j  (x)F{t) 

'  o'  9,,  (♦)  -  - -  (15) 

where  M  is  the  generalized  moss  and  u,,  (  27tf,,  )  the  notural  frequency  of  toe 
I'-th  mode.  The  steady-state  solution  of  Equotion  (15)  when  F(t)  is  harmonic  of 
frequency  f  (  Ui/Tr,)  is 
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4),,  (x)  F(t)/  M., 


q..  (t) 


u  u,, 


(x)  H„(f)  F(t) 

t' 


(16) 


where  k,,  (  u),,  M^,)  is  the  generalized  stiffness  of  tfie  i'  -th  mode  and  where 


H,  (f) 


1 


1  ^ 

o  U. 


To  obtain  the  modal  acceleration  response,  replace  q,,  (*)  by  -  q^,  (t)/u  in  Equation 
(16)  giving 


r  2 


n,0)  '1^-^—  <f>.  (x)  H,,(f)  •  r(t) 

k 


(17) 


The  total  acceleration  V,'ft)  of  tlie  beam  at  y  can  be  found  by  substi tijting  Equotiori 
(17)  into  Equation  (I'D  giving 


W(f)  a  (f)  F(t) 


where  the  beam  transfer  function  a  (f)  is 


2  11.  (D 

a  (f)  ~(2Trf)  /  j  — ti,  (y1 


(18) 


Ttrus  it  is  seen  tliat  tire  force-to-acceleration  beam  transfer  functic<n  can  lie  expressed 
us  on  infinite  series  of  normal  rriode  functions  for  the  beam. 
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Mcon  Square  Acceleration  Response 


Returning  now  to  tlie  general  expression  for  the  net  acceleration  response  W^(t)  In 

Equation  (4),  It  is  instructive  to  write  this  equation  in  terms  of  flie  absolute  values 

of  a  (f,  )  and  F  ,  (t),  the  phase  angles  6  (f,  )  introduced  by  the  structural  trans- 
nr  k  nk  nr  k 

mission  paths,  and  the  elative  phase  angles  of  the  applied  forces.  The  applied 
force  component  F  |^(t)  is  harmonic  of  frequency  fj^  and  con  therefore  be  written 
in  the  complex  form 


(t) 


'•^nk  ^'^k* 


(19) 


where  jf  denotes  the  amplitude  and  ^  ^  the  phase  of  Substituting 

Equation  (7)  for  a  (fj^)  and  Equation  (19)  for  1^^|^(*)  Equation  (4)  gives  the 

•  « 

following  expression  for  W  (t) 


W^(t) 


E  E 


2Trif,  t 
k 

p 


Equation  (20)  shows  that  the  net  occeleratlon  depends  upon  both  the  phase  angle 
chonges  which  take  place  along  the  structurol  tronsmisslon  paths  and  the  relotive 
phase  angifs  of  tfre  applied  forces.  When  the  mean-square  value  (t)  of  the 
net  acceleration  is  to  be  determined,  It  Is  necessary  to  take  the  time  average  of  the 
squore  of  the  real  value  of  W^(t),  that  is 


W^t)  - 


<!'  (21) 


whe 


re 


(f  {*  ('>j 


denotes  the  real  part  of  W  (t) 
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Substituting  Equation  (?0)  into  Equation  (21),  and  noting  that  response  components 
of  diffeient  frequencies  are  uncorrelated,  gives  the  following  expression  for  the 
mean-square  value  of  W  (t). 


WMt) 


2  I 


fn  /  m) 


*  cos  0  (f  )  -  9  (f,  )  * 

nr  k  mt  k 


Equation  (22)  shows  that  the  mean-square  volue  of  the  net  occelerotion  response  is 
equal  to  the  summation  of  mean-square  response  components  and  crcsscorrelotions 
of  response  comptvonents .  Tlie  quantity 


1  (f  \ 

T  ('i  ) 

/  nr  k 


is  the  mean-squore  response  at  "y  due  to  on 


inpi/t  force  at  x  whicli  is  liarmonic  at  frequency  f,  and  whicli  has  an  omplitude  of 
n  '  k 

jr  .  The  terms  within  the  second  set  of  summations  in  Equation  (22)  are  the 

crosscorrelations  of  responses  at  y  ,  for  the  same  frequency  f^ ,  but  for  excitations 

at  points  "x*  and  x 
n  m 


It  is  seen  tliaf  tlie  magnitudes  of  the  crosscorrelotion  terms  are  dependent  upon  the 

r  1 

pilose  angle  differences  0  (f,  )  -  0  (f,  )  |  which  occur  for  different  structural 

^nrk  mrkj|- 

paths  and  upon  the  phase  ongle  differences  rT*  p  “  ‘1’  p  I  between  the  forces 

applied  at  x  and  x  .  Clearly  the  crosscorrelation  terms  in  Equation  (22)  are 
n  m 

exactly  equal  to  zero  only  for  the  two  special  cases  in  which  either  there  is  only 
one  applied  force  or  tlm  phase  angles  ate  such  that  the  argument  of  the  cosine  term 
in  Equation  (?'?)  is  equal  to  r;/ 2,  or  an  odd  multiple  thereof.  In  the  sp<*cial  case 


22 


n 


when  oil  of  the  forces  ore  in  phase  end  all  transmission  paths  cause  the  same  phase 
angle,  the  argument  of  the  cosine  in  Equotion  (22)  is  equol  to  zero,  ond  *he  mean- 
square  response  expression  reduces  to  the  form 


i  “ 

• 

r  i  X 

L 

la  (U 

I  nr  k 

■  rnkl 

!  h 

n 

1  1 

n  2 


(23) 


Further,  it  is  to  be  noted  in  Equation  (22)  that  the  mean-square  terms  are  positive  and 

thus  each  mean-squate  response  component  adds  to  the  value  of  W ^t ) .  The  presence 

of  the  cosine  foctor  in  the  crosscorrelation  terms  can  lead  to  negative  valued  cross- 

correlotlon  terms,  so  that  In  some  instances,  crosscorrelation  terms  tend  to  reduce 

the  magnitude  of  (t) .  In  fact,  it  is  possible  for  the  ctosscorrrelotion  terms  to  be 

of  such  a  magnitude  and  sign  that  the  total  mean-square  response  is  equal  to  zero. 

This  can  be  demonstrated  for  the  following  elementary  example.  Consider  the  rigid 

bar,  shown  in  Figure  4,  which  has  a  large  number  of  applied  forces  F  (t),  all  of 

n 

which  are  harmonic  at  the  some  frequency,  and  have  the  some  amplitude  jFj  . 

The  bar  is  assumed  to  be  free  to  translate  but  constiained  in  pitch.  The  transfer 
functions  a  (f.  )  for  this  case  are  equol  to  1/M  where  M  is  the  mass  of  the  beom; 
the  structural  phase  angles  he  set  equol  to  zero;  and  the  force  amplitudes 

F^j^j  can  be  replaced  by  |f  [  .  Thus  if  there  ore  N  such  forces,  the  mean-square 
acceleration  response  V/^  (t)  in  Equation  (12)  can  be  written  in  the  form 


(t) 


2M 


N  4 


N  N 

E.  E 


m 


cos  (<;>-({)) 


(24) 


(n  /  m) 

where  (j)  and  rf  denote  phase  angles  of  the  applied  forces.  If  the  phase  angles  * 
hove  a  uniform  probability  of  occurrence  In  the  range  from  0  to  2n ,  then  the  double 
summation  in  Equotion  (24)  is  approximotely  equal  to  -N;  and  it  is  easily  shown  that 
when  the  phase  angles  are  uniformly  distributed  between  0  and  2it  ,  this  double 
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summation  is  exactly  equal  to  -N  when  N  is  an  even  integer .  Thus,  in  this  case 
the  summation  of  all  of  the  crosscorrelation  terms  is  equal  in  magnitude  to  the 
summation  of  oil  of  the  mean-square  terms  in  Equation  (22)  end  there  is  no  net 
lesponse  . 

Cros^oi  re  lot  ion  of  Acceleration  Response 

In  anolyring  vibration  responses  of  structural  systems,  it  Is  often  of  interest  to 

determine  the  crosscorrelation  of  tlie  responses  at  two  different  points  of  the 

structure  for  o  given  loading.  Consider  once  again  that  point  forces  Of® 

applied  to  the  structure  at  "x*  ,  and  consider  the  net  acceleration  responses 

^  n 

W  (t)  at  "y*  and  W  (t)  at  "y*  ■  expression  for  W  (t)  is  given  l>y  Isrjualion 
r  r  s  s  ^  f 

(4)  and  the  comporable  equation  for  W^(t)  is 

K  m 

The  crosscorrelation  of  these  two  acceleration  responses  is  denoted  by  rj(f)  ond 
can  be  defined  os  follows 


1 


dt 


(26) 


wlieie 


•  » 

denotes  the  real  part.  It  is  clear  from  the  form  of  W  (t)^in  Equal.' 


on 


(20)  that  the  reol  part  of  W  (t)  is 


r?[w(.)]  E  E  kA: 

K  n  I 


''nkj  '  0„<'k’  *  ♦nkj 


(27) 


Similarly  for  the  real  part  of  W  (f) 
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(28) 


fmk[  [^'k' " '  ®ns"k>  '  ♦„!< 


SuhstI futing  Equations  {?7)  and  (28)  into  Equation  (26),  peifoiming  the  necessary 
integrations,  and  noting  tliat  the  response  components  of  different  frequencies  are 
uncorrelated,  gives  the  following  expression  for  the  crosscorrelation  function  I’  (r) 


T 


ESS 


n  m 


a  (f,  )  •  a  (f.)  •  h  , 
nr  k  ms  k  nk 


mk 


(29) 


cos 


^"'k’  '«J'k'  -  8.,<'k>  '  t„k  -»nk 


] 


It  Is  seen  that  Equation  (29)  reduces  to  the  meon-squate  response  given  by  Equation 

(22)  when  "y*  ,  "y*  and  ^  ■  0. 

s 


POINT  FORCE  RAN  DOM  EXCITATION 

Definitions  and  Approach 

Tire  acceleration  trs|)onse  W^(t)  at  the  point  on  a  structure  subjected  to 
randomly  varying  forces  f  ft)  applied  at  the  points  con  be  described  by  the 
methods  of  hormonic  onrtlysis. 

Consider  tire  excitation  forces  F^(t)  and  the  response  W^(t)  to  be  stationary 
random  processes,  and  assume  that  they  began  at  t  ‘  -co .  In  a  strict  sense,  the 
Fourier  transform  of  W^(t)  rcjnnot  be  defined  over  an  rnfinite  time  Interval  for  a 
stationary  process .  However ,  the  Fourier  transform  of  W^(t)  can  be  defined  for 
that  part  of  W^(t)  over  the  interval  -T/2  <  t  <  T/2,  assirming,  of  course,  fjrat 
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w 


V  (t)  is  zero  for  all  other  times,  hence,  letting  W  (t)  represent  the  response 
t  ■  2 

in  the  above  time  interval,  the  mean-square  value  W  of  W^.|.(t) 


iut  transform 

« 

T/2 

f 

-fr 

M/7 

oo 

/ 

r 

-oo 

1 

oo 

/ 

*.T< 

•00 

oo 

/ 

W,(i 

-oo 

OD 

I 

-CO 


A_(f)  e'^^^dfdt 

tl 


Integrating  first  with  respect  to 

time  gives 

•2  1 

V/.(f)  -r 

iT  ' 

7 

.  i?iTft 

-00 

-CD 

1 

T 

oo 

/ 

A*  (f)  df 

1 1 

-oo 

where  A  *(f)  is  the  compi 
rT 

ex  conjugate  of  A  ^(f),  so  that 

Wf,(0  f 

oo 

/■ 

A  ^(f) 

f  1 

^  d. 

Let  T  ♦•'TO,  cinrl  an  expies'.ion  for  llm  mean  squore  value  of  \V  (f)  results 
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w^^(0 


7 


(f)  df 


(30) 


•  • 

where  Sj.(f)  is  defined  as  the  spectral  density  (S  .D  .)  of  the  response  W^(t) . 

Trr~  /  V  .  . 

Therefore,  in  computing  tire  mean-square  value  W  (t)  of  the  response  W^(t)  it  is 

seen  that  the  spectral  density  S^(0  must  he  determined.  The  spectral  density 

S  (f)  indicates  the  distribution  of  the  hoimonic  content  of  W  (t). 
r  f 


The  S  .D  .  S  (f),  howevei ,  can  be  determined  from  a  knowledge  of  the  autocorrelation 
r 

function  !'(■»)  of  response  W  (t)  and  the  use  of  the  Wiener  -  Khincbln 

relotlons.  The  autocorrelation  function  )  is  defined  in  terms  of  tire  response 

W  (t)  as 
r 


W  (t)  W  (t  t  T  )  dt. 
r  r 


(31) 


In  general,  the  notation  for  time  overoge  will  be  simplified  as  follows 


l-(-r  )  VV  (t)  W  (t  '1  )  .  (32) 

r  r  r 

The  autocorrelation  1'  (r)  and  the  spectral  density  5^(f)  of  the  response  W  (t) 
are  Fourier  transform  fxrirs,  so  that 
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s  (0  2 


7  ■' 


,  >  -i2TlfT 

(■r  )  e  di  , 


1  ^ 

'  >  T  /  ^ 


(f)  d(. 


and  fhe^e  exnressions  are  known  as  the  Vviener-Khinchin  relations. 

In  Older  to  develop  further  these  expressions  for  response,  it  is  necessary  to  relate 

the  response  to  the  excitation  forces  using  the  tronsfer  characteristics  of  the  structure, 

The  force-to-accelerotlon  transfer  function  n  (f)  and  the  unit  impulse  restxsnse 

nr 

function  h  (t)  completely  describe  the  dynamical  properties  of  the  structural  paths 

between  the  points  x*  ,  wh<'re  the  forces  F  (t)  Oi  e  applied,  and  the  points 

n  n 

_ ^  •  • 

,  wtiere  the  responses  are  lieing  considered.  Tire  response  W  (t) 

can  be  defined  in  terms  of  the  unit-impulse  response  function  h  (t)  end  the  opplied 
forces  F  (t)  by  use  of  the  convolution  integral 


W  (t) 


h  (t)  F  (t--r)dT 
nr  n 


Further,  tfie  transfer  function  a  (f)  is  defined  in  terms  of  the  Fourier  transform  of 

nr 

the  unit-impulse  response  functiori  li  (t)  , 

nr 


a  (f)  -(2nf7  f  h  (t  d-r  . 

nr  #  nr 

*t) 
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Substituting  the  inverse  Fourier  transform 


7/n 

a  (f)e 

-  df 


(36) 


info  equation  (34)^  flip  expression  for  the  acceleration  response  becomes 


a  (f)  e 
nr 

(2rTf)^ 


F(t  -  T  )  df  dr  • 


(37 


However,  a  more  convenient  opproach  wculd  involve  the  development  of  an  ex¬ 
pression  foi  the  arjtocorrelation  f’  (t  )  of  W^(t)  in  terms  of  the  autocorrelations 

r  (t  ),  r  (t  ),  tlie  crosscorrelofions  j'  (i  )  of  the  applied  forces  F  (t),  F  (t) 
m  n  mn  m  n 

end  the  unit-impulse  transfer  functions  h  (t)  and  h  (t).  Tlie  crosscoi relation 

ml  nr 

I  (t  )  of  the  forces  F  ft)  and  F  (f)  is  defined  as 
mn  m  n 


r  (t  )  F  (t)  F  ft  It  )  .  (33) 

mn  m  n 

The  spectral  density  S  (f)  of  W  (t)  con  then  be  related  to  the  spectial  densities 

r  r 

S  (f),  S  (f),  the  cross-spectral  densities  S  ff)  of  the  forces  F  (t),  F  (t),  and 
m  n  mn  m  n 

the  transfer  functions  a  (f),  a  (f)  and  their  complex  conjugates.  The  relation 

mr  ni 

Isetween  crosscorrelotion  and  C.S.D.  is 


s 

mn 


(f) 


,  .  -l2nfT  , 
r  (t  )  e  dT 

mn 


(39) 
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In  anolyzing  the  vibrations  of  complex  structure  it  is  often  of  Interest  to  determine 

the  cfrosscorrelotlon  T  (i  )  and  cross -spectral  densities  S  (f)  of  the  acceleration 

rs  rs 

responses  and  V)^(f)  at  the  two  points  and  .  These  crosscorrelctlons 

are  defined  as 


I'  (7  )  W  (t)  W  (t  ^  T  )  .  (40) 

rs  f  s 

artd  the  C  .S  .D  .  of  the  two  responses  is  defined  as 


CD 

S  (f)  '  2  /  p  (7  )  d7  (41) 

rs  rs 

For  example,  knowledge  of  the  above  relations  allows  a  complete  description  of  the 

responses  at  various  points  witliln  a  structural  compartment  of  a  vehicle  excited  by  multiple 

external  forces.  Also,  a  particuior  structural  path  moy  be  traced  by  use  of  cross- 

correlations  1'  (r)  or  cross -spec'.tal  densities  S,  .(f)  of  the  difference  of 

(r  -  s)  ( r  -  s; 

the  responses  Vv'  (t)  and  W  (t)  at  two  different  points  ‘y  and  y  .  Such  cross¬ 
correlations  and  cross-spectral  densities  are  defined  by  the  equations 


'(r  -s) 


(r) 


W  (t)  -  W  (t)  W  (t  '1  )  -  w  (t  <7  ) 
r  sir  s  . 


(42) 


S-  /f)  2 

(r  -  s) 


'  (r  -  ' 


-i  ?Ttf7 


d7 


(43) 


The  various  expressions  given  above  are  developed  In  detail  in  the  discussions  to 
follow.  , 
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Response  at  o  Single  Point  Due  to  o  Single  Loading 

A .  Autocorrelation  ond  Spcctrol  Density  of  Excitotlon  and  Re'.ponse 

•  • 

The  autocorrelation  I'  ^  t  )  of  the  acceleration  response  W^(t)  at  a  single  point 
due  to  a  single  force  F  (t)  at  a  point  x  is  obtained  by  substituting  Equation 
(34)  into  Equation  (32)  and  rearranging 


r  (7  )  -  (t)  W  (t  1 1  ) 

r  r  r 


(?rTf)' 


h  (7.) 
nr  1 


00 


(2nf)^h  (tJ 
nr  / 


F  (t  - 

n 


)Fn(. 


For  stationary  processes,  n  change  in  the  origin  of  t  will  make  nc  char>ge  in  the  meaned 
quantity .  Therefore 


1(7  ) 

r 


(?nf) 


7h  (. )  J 


h  (i„) 
nr  / 


■  (7, 

n  1 


(44) 


This  relationship  between  t!ie  autocorrelation  p  (t  )  of  the  response  ond  the  auto¬ 
correlation  I'  (t  )  of  the  force  is  not  a  simple  one  .  However,  the  difficult 
n 

relotionshia  above  will  result  In  a  simple  relation  between  the  spectral  density 
S  (f)  of  the  response  and  the  'oectrol  density  S  (f)  of  the  input  force  by  using 
the  Fourier  transform  relarion  of  Equation  (33), 


no 

2 


5(f)  2  /  L(T)e"'^^^dT 


-OO 


00  00  00 

^  2  2  2  '  n<’) 


2^'1 


-i2itfT  , 
e  d 


31 


Rearranging  »iiis  integral  give^ 

S  (,) 


(2nf)‘h  dT,  /  (2«.f)^h  (t  ) 

0  ’  •'O  ^  ^ 


7  '.v 

-CO 


-i2nf(T  -  1-  ♦  T  ) 

‘-r  )  e  d(T^  -  ♦■r  )  , 


and,  finally,  from  Equations  (35)  and  (33),  this  expression  reduces  to 


or 


S(f)-  a*  (0  o„(f)S  (f) 

r  nr  nr  n 


S^(f) 


S  (f) 

n 


(45) 


B.  Response  of  Single  Degree-of-Freedom  System 

Consider  now  the  single  degree  of  freedom  system  shown  in  Figure  2  end  assume 
that  the  applied  force  is  random.  In  order  to  maintain  the  same  nototion  as  thnt 
used  above,  the  spectral  density  of  the  acceleration  response  of  the  mass  is  denoted 
by  W(f)  and  the  spectral  density  of  the  applied  force  Is  denoted  by  S^(f) .  Sub¬ 
stituting  Equation  (9)  Into  Equation  (45)  gives  the  following  relotionship  between 
the  spectral  densities  of  response  ond  excitation 


S  ff) 


(2nf)'’  S  (H/ls^ 
n 


(46) 


where  f  is  the  resonant  frequency  and  k  the  spring  constant  for  the  single 
o 

degree-of-freedom  system. 
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The  mean-square  value  of  the  responi'  W  (t)  can  be  obtained  by  Integrating 
Equation  (46)  over  all  frequencies, 

_ _  00 

(t)  -  f  S^(f)  df  •  (47) 

0 

Now,  if  the  excitation  force  F  (t)  is  white  noise,  then  the  spectral  density 

n 

S^(f)  can  be  replaced  by  the  constant  S;  ond  for  this  case,  Equation  (46) 


S  (f) 

r 


4  2 

(2irf)  S/k"^ 


2 

T  / 

r( 

“o). 

b) 

(49) 


And,  hence,  the  mean-square  acceleratlor,  response  for  the  single  degree-of-freedom 
system  Is  given  by  the  equation, 


which  upon  using  residue  theory  gives 


(49) 


C.  Crosscorteiation  ond  Cross-Spectrol  Density  of  Excitation  orxt  Response 
The  crosscorteiation  J'  (t  )  of  the  single  Input  force  F  (t)  and  the  response 

•t  «  ^ 

W  (t)  can  be  obtoined  by  substituting  Equation  (34) 


••  2 

wMt) 

r 


SO 


2  3 

m  (?nn 
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I-  (t)  =  ‘f  (t)  VJ  (t  tT  ) 
nr  nr 


■^(t)  /‘-{2rrf)  h^^(T,)  F^(t+T-r,)  dt, 

^r\ 


'0 


~(2t.f)^ 


f  ’-■'l'  '*’1 

*^t\ 


r  (t  ) 

nr 


7 


(2nf)  I  h  (t.)  ['  (t  -  T.)  d-T, 

*  nr  I  n  '  1'  1 


(50) 


which  gives  the  crosscorrelation  )  in  terms  of  the  unit-impulse  response 

function  h  (t)  ond  the  outocorrelation  1'  (i  )  of  the  input  force  F  (t), 
nr  n  n 


The  cross-spectral  density  S^^(f)  of  the  input  force  F^(t)  and  the  response  W^(t) 
is  defined,  as  in  Equation  (39),  as 


00 

"  2  /  r  (-r  )e"'^^"  dT 

nr  j  nr 

-oo 

and,  from  Equation  (50) 


34 


00^  r  00 

7  7 

*40  [  *'0 


'  (2nf)‘ 


-i2lTfT^ 


7  'n'’  -  ’V 


■i2nf(T  -  T 


d(T  -  T^) 


S  (f)  -  a(f)  S  (f) 

nr  nr  n 


The  cross -spectral  density  S  ^(f)  of  the  input  force  ^^^(0  ond  the  response  'W^{t) 
is  simply  the  product  of  the  transfer  function  a  (f)  and  the  spectral  density 


5  (f)  of  the  Input  force  F  (t). 
n  n 


The  ctoss-spectrol  density  S  (f)  in  Equation  (51)  can  olso  be  written  in  terms  of 

'  nr  ^ 

the  spectral  density  S^(f)  of  the  response  W^(t).  The  complex  conjugate  of  Equation 
(51)  is 


S‘(f)  a*  (f)  5  (f) 
nr  nr  n 


since  5^(f)  is  a  real  quantity.  Multiplying  Equation  (51)  by  a  *  (f) 


a  •  (f)  S  (f)  -  a*  (f)  a  (f)  S  (f) 
nr  nr  nr  nr  n 


'l'-nr<'T  'n<'> 


-  S(f) 
r 
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Hence, 


S 

nt 


(f)  •- 


1 

a-  (f) 

nr 


S  (f) 


r 


(53) 


The  spectral  density  of  white  noise  is  o  constant  S.  Thus,  the  cross-spectral  density 
of  the  acceleration  response  V/  (t)  of  a  point  on  a  structure  excited  by  white 
noise  is  from  Equation  (51) 

S  (f)  -  a  (f)  S  (54) 

nr  nr 

The  cross-spectral  density  S  (f)  has  the  same  form  as  the  tronsfer  function  a  (f)- 


Tire  crosscorrelation  function  I'  (t  )  of  tire  white  noise  excitation  and  the  response 
W^(t),  using  tire  definitlotr  of  the  autocorrelation  I'  (t  )  of  white  noise  given  by 
Ecpjation  (A.  10)  Appendix  A  ,  os 


I-  (-r  )  -1  S6(.  ), 

n  / 

is  written  as 

7 

I  (-r  )  (2nf)^  /  h  (tJ  !'  (-.  -  -r.)  dT, 

nr  nr  1  n  11 

2  ^ 

■  /  h  (t  )  5  6(-r  -  1.)  di 

2  nr  1  11 

I-  (i  )  2..^f^  S  h  ).  (55) 

nr  nr 
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For  white  noise  excitation,  therefore,  the  crosscorrelation  of  excitation  and 
response  has  the  form  of  the  unit  impulse  response  function  )• 

The  above  result  has  opplicotion  os  o  good  experimental  method  for  determining 
the  vibration  transfer  characteristics  of  structures. 


Response  of  Multiple  Points  Due  to  Multiple  Loadings 

A.  Autocorrelation  and  Spectral  Density  of  the  Response  to  Multiple  Loadings 

•  » 

Tlie  autocorrelation  )  of  the  acceleration  response  W^(r)  at  a  single  point 
on  a  structure  subjected  to  three  forces  F^(t),  and  F^Ct)  v/hich  ore  applied 

at  points  ,  'x*^  and"!^,  respectively,  is  expended  from  Equation  (32)  by  sub¬ 
stituting  from  Equation  (34) 


r(T  )  -  W  (f)'//  (t  -  T  ) 

r  r  r 

(2,0'’  7k,(’|)  F,('  -  •>,)  '  Fjd  FjO  d, 

n)  J 

— -  _ 

■  /  h<'  -’2*  ■'  >  '  '’2,<’2>''2"  -'2'  'S,'’2’f3<' 


and  multiplying  and  rearranging. 


3  3 


r(0  (2nf) 


EE? 

m  -  1  n  -  1  / 


\XJ 

f  f  (0  F  (t*  -r  -1-  +T  )  dT 

'  /  nr  /  m  n  I  2 
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And,  since 


r  (t,  -  T,  *  T  )  ■  F  (0  F  (t  ‘T,  -  ^  T  ) 

mn  1  2  m  n  1  2 


t!)cn , 


(0  (2lTf) 


3 _  3  oo^  oo, 

'*  Z]  /  J  '’2'  '^’2  ‘‘’I 

m  1  n  1  '  0  0 


The  autocorielation  function  I'  (t  )  of  the  response  to  three  input  random  forces  ’s 

now  written  in  teims  of  the  three  outocorielotion  functions  I'  (t  )  aid  the  six 

n 

ciosscotteintion  fi/nctions  I'  (i  )  of  t!ie  applied  forces  F  (t)  and  F  (t)  and  the 

mn  m  n 

appiopiiote  unit-impulse  lespon'.e  functions  h  {^  ). 

mi 


The  spectial  density  S^(f)  of  the  lesponse  VJ  (t)  is  deteimined  from  the  Fourier 
tinnsfrjim  i-lntion 


S  ff)  2 

r 


r^{-r  )  di 


CD  GD^  OU 

f  /  (■r,-T_tT)di  d7,  e'^^^di 

/  /  mr  I  f  nr  2  mn  1  2  2  1 

^-OD  ''n  *^0 


/  h  dT.  fu  iije 

/  mi  1  1  /  nr  2  2 


-i  ?nf  ~  "^2  '  ^  ^ 

-•'2'  d(,, 
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3  3 

S  (0  -  (0  °  (0  S  (f), 

r  Amc/,  /arcstf,  mr  nr  mn 

m  -  I  .  n  ■  I 


where 


S  (f) 
mn 


/  r  (t 


.  -i2nfr  , 
)e  OT 


is  the  cross-spectral  density  of  the  applied  forces  F  (t)  and  F  (t). 

/  rr  m  n 


Now,  In  general  the  autocorrelation  T  (t  )  and  spectral  density  S  (f)  of  W  (t)  to 

n  r  r 


N  applied  forces  ore,  respectively, 


J  F'  (t,  -  ^  t  )  dT_  dT,  (3 

2  mn  I  2  2  1 


N  N 

m  l  n  I 


To  illustrate  the  effect  of  crosscorrelation  on  the  response  of  a  point  In  a  structure, 
consider  the  previous  case  of  the  three  forces,  where  (1)  the  three  forces  are  indepeiide 
and  no  crosscorrelotions  ex'st  between  the  forces,  and  (2)  the  three  forces  are  perfectly 
correlated  and  have  identical  spectral  densities.  In  the  first  case,  the  spectral  density 
the  response  is  - 


S^f)  = 


Z. 


a  (f)  S  (f)  . 
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Now,  let  5.(f)  S_(f)  ^  S„(f)  ■  S,  and  the  spectral  density  of  the  response  in 

I  C 

the  second  case  becomes 


where  0  are  the  relative  phase  angles  between  each  pair  of  transfer  functions  ot  the 

mil 

patficulor  frequency  f. 

Therefore,  It  Is  seen  that  when  there  Is  correlation  between  the  opplled  forces  the 
spectral  density  of  the  response  at  any  frequency  depends  upon  the  relative  phase 
angles  between  the  transfer  firnctlons  of  the  structural  paths.  And,  wlien  no 
correlation  exists,  the  relative  phase  angles  have  no  effect  on  the  response  magni tude  . 

B.  Crosscorrelation  and  Cross-Spectral  Density  of  Two  Responses  to  Multiple  Loadings 

Consider  the  responses  W  (t)  and  W  (t)  at  points  and  to  three  input  random 

r  s  r  s 

forces  F^{f),  F^ft)  and  F2(t)  applied  at  points  respectively.  The 

crosscorrelofion  function  I'  (t  )  of  the  responses  W  (t)  and  W  (t)  can  be  obtained 

rs  r  s 

by  expanding  the  responses  ns  before  from  Equotion  (34) 
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r  (t )  -  w  (0  w  ( t  ) 

rs  r  s 


•  7  h|^(,2)F,(l-.j  M  )<  thj^hjjFjO  -Tj 

*^f\ 


^T  )  dT, 


multi  pi  yin-3  rearranging,  gives 


3  3 


1  (t)  -  (2^»0 
rs  ' 


1  n - 1 


OD  03, 

fh  (t.)  /  h  (tJ  r  (t,  -  T  t  T  )  dT  dr,  (62) 

I  mr  I  I  ns  /  mn  \  Z  Z  \ 

^(\ 


Similarly,  tire  cross -spectral  density  S  ^{f,'  is 


3  3 


rs  '  mr  ns  mn 

m  Ini 


Obviously  the  extension  of  Equations  (62)  and  (63)  to  the  cose  of  K  random  loadings 


N  N 


I  (r)  -  (2nf) 


b  (t,) 


1  n  1 


oo 

f 


>7)  I'  (f  1  ~  ■'n  '  T)dT_  dT, 

2  mn  1  2  '21 


N  N 


S  (f)  7  ^  a  (f)  a  (f)  S  (f), 

rs  iLi^,  rnr  ns  mn 

mini 
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This  completes  the  general  theory  describing  the  response  of  multiple  structural 

points  excited  by  multiple  rondom  loadings  applied  at  other  structural  points.  A 

description  of  the  responses  at  many  internol  points  of  a  st.uctural  segment,  such  as, 

a  compartment,  etc.,  relating  each  response  VV^(t)  to  each  response  'W^(t),  and  to 

each  input  force  F  (t)  as  well,  can  be  obtained, 
m 

C.  Crosscorrelation  and  Cioss -Spectrol  Density  of  Beom  Responses 

Consider  the  uniform  beam  of  length  L  ,  sliown  in  Figure  5  witli  normal  modes 
(j)  (y)  and  natural  frequencies  f  .  The  generalized  mass  In  each  mode  Is 


/ » (7>  rn  d^ 


where  m  is  the  mass  per  unit  length  of  the  beam.  The  cross-spectral  density  of 
the  acceleration  responses  W  (t)  and  W  (t)  at  and~^  due  to  an  input  random 


force  F  (t)  at  y  is, 
n  n 


S  (f)  a  (f)  a  (f)  S  (f). 
rs  nr  ns  n 


Writing  the  transfer  function  in  complex  form 


a  (f) 
mt 


z 


<f>  (7  )  t  (y  ) 
m  q  ' 


r?  .  r 

f  -ft  icf 

^  2rT  M 


Et  (^  }  <p  (y  ) 

^q  'm  q  r 


■(f  '  (f) 

\q/  q  q  ' 
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wlilfl;  i  ill  ()PM(>inl  i-oriplc/  fincl  cfimKit  t  '■  si  mpl  i  fl  i^rl  nily  bfcauso  of  llif'  orodiyct 

|p|  rns  . 


Sirnplificoti  'MS  i<''-'/lt,  tiov.rvci,  fio'i'  tli*'  nsS'i'njsl  ions  tlicil  tlio  stiuctural  rioniplnq  is 
snioll  one)  dint  llin  i csfvinst'-  D«'“nl-s  om  v,r>ll  scpciiotecj.  Tlmn  tho  piodnct  torms  wlioio 
q  j  piodominoto  ovrr  tfiose  v-lioio  q  /  j  /  <111(1,  tlu'rofort',  llu’  l.'ilttM'  can  lie  negl('cto 
Honco , 
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D.  Response  at  Twj  Points  Due  to  Two  Inpot  Fo  C  -s 

Now,  consiHer  the  piohler.i  of  computing  the  crosscorrelation  of  two  acceleration 

responses  W  (t)  ond  W^(t)  at  points  and  on  a  structure  '  -cited  by  two 

input  forces  F  (t)  and  F  (t)  applied  at  points  x  and"?  .  The  uni t-impulse 

response  functions  h  (t)  ,  h  (t),  h  (t)  and  h  (t)  describe  the  dynamic  character- 
mr  ms  nr  ns 


Istics  of  the  structure  betv/een  the  points  y  ,  y  and  x  ,  x 

r  s  m  n 

schemotically  illustrates  the  problem. 


Figure  6 


From  Mriiiation  (fil)  the  form  of  tlw'  cro.s.‘roorrolation  funution 


r  (r)  -  (2rrf) 
rs 


2  2  '' 


1  /  ns  2  mn  1  2 


1  n  1 


Hence , 


i;^(r)  (2^f) 


/  b  (rj  /  h  (tJ  I  (i.  -  T,  t  T  )  dT- dT, 
J  tr\r  \  J  ms  2  mn  1  2  2  1 

-oo 


+  /  *’  (■’i)  /  ^  ''  i  ''')  '  '  '  ■*■^1 

#  mr  I  f  ns  2  mn  12  2  1 


no  OD 

^  /  b  (t.)  C  b  (t 

I  nr  1  #  ms  2 


)  F  (t,  -  T.,  ♦  T  )  di-  d-r, 
nm  12  21 


*  j  '’nr<’l>  y  '’2”’  (73) 


Now,  let  7^  *  '  ^1  EquatIon(73) 
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The  crosscoitelation 


I  (  T  )  of  the  two  responses  is,  from  Equation  (75). 
nnr,ms 


r 

mr ,  ms 


(t) 


h  (t)  h  (t  '  T  )  dt 
mr  ms 


(79) 


If  T  -^0,  then  from  Equations  ^77)  and  (78)  the  abose  integrand  equals  zero  in  the 
range  -cn  <  t  <  0  so  that  Equation  (79)  reduces  to 


I' 

mr 


(t  ) 


rns 


h  (t)h  (t+T)dt,  t^O 

mr  ms 


(80) 


When  T  'I  0,  then  Equation  '79)  is  wiitlen 


I  (-hP 

‘mr,ms  ' 


00^ 

/  h  (t)  h  (t  -  1  7  I  )  dt,  t  <  0 

j  mr  ms 
‘-00 


(31) 


where,  from  Equotions  (77)  end  (78)  it  is  seen  that  the  integrand  above  is  zero  in 
the  rancie  -oo  \  "  |i  (  so  mat  Equation  (31)  reduces  to 


I  (-hP 

mt,ms  '  ' 


h  (t)h  (t-|7|)dt,  74  0 

mr  ms 


And,  leOing 


I  7  I  ,  gives 
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I  (-■  hi  )  /  h  (t  '  |i|  )  fi 

rni  ,  ms  '  #  mr  ms 

‘0 


h  (t  *  |i|  )  fi  (0  dt  ,  T  ^  C. 


Equations  (80)  and  (82)  orp  used  to  compute  the  crosscorrolation  (t  )  of  the 

mr ,  ms 

responses  h  (f),  h  (t)  of  two  slnnle  deqrpe  of  freedom  systerris  to  tfie  Impulse 
mr  ms  ' 

pycitotlon  F 

rn 


Substituting  Eqijatlons  (77)  and  (73)  Into  Equation  (30) 


.  -CTrrfT 

F  e  ^  ^  oo  -?n(s  f  tF  f  ) 

-  >  m  /  t  r  s  s 

1  n  )  - s -  /  e 

M  m  on  ff  (7  J 

r  $  r  $  0 


sin  2tt  ft  t  sin  ?TT  <T^{t  *  t  )  dt,  t  0 


and  expending  by  using  e  tt igonometi I c  identity. 


-  -  f  ?n  f  T 

5  S 

F  e 

I  (t  )  - 5 - 

r  s  r  s 


cos  2Tr 


"  ^  J 


00  -?TT  (  f  f  f  c  f  ) 
/"  r  r  s  s 


sin  ?n  a  t  sin  2it  a  t  dt 
r  s 


CD  -2it(  n^) 

‘  sin  ?n  ,T  T  /  ft  sin  ?n  (T  t  cos  2tt  <r  t  dt  ,  t  '  0 

i  /  r  s' 
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Integ  atlrig  ntid  collec  tiricj  tot'in 


I’  (t 

mr ,  ms 


(  ?nf  T  r 
s  5 


a  cos  /n,TT  •hsin?nr77 
s  s 


0  (‘^3) 


whpi  ( 


2ii  Y  F 


.2 


32*t  t\\  IA  t  r 

I  SIS 


(a  -  a  Y 

t  s 


♦  (o  'a 


r  s 


32ii  M  K\  n  n 
I  s  f  s 


a  -  a 
I  s 


a  '  a 
r  s 


‘  (rj  -  a  ) 
'  I  s 


2  .  .2 
'  (a  >  0  ) 
r  s 


>  ,  f  •  I  f 

II  s  s 


Eqi.inti^'M  F--,2)  sIhts'-s  lii'it  t!if  c..''i r'*5p''ndMi(!  I'yprf’ssion  foi  nf.-qatlvc  values  of  i  con  he 
sshtaiivv!  Isy  inteirlirMVjin')  »'ri«-  and  in  Equation  (33^  OrKj  ifplcicinq  t  witldrl 
The  oh'v  io'iS  lesul  t  is 


(-Ml 

-  C  2tt  f  111 

r  r  * 

n 

c  cos  ?n  a  ItI 

'  d  s  i  11  2n  a  I  T I 

I  ,  T  ' 

rri.fTis  1  ' 

t  '  ' 

r  '  '  J 

vv  n  *•'  I  e 

0 

2^Y 

'  m 

1 

1 

3?it  M  M  a  a 

I  SIS 

2  2 
y  ‘  (o  -  0  ) 

L  '  SI 

"■  2  2  ■ 
y  *  (0  ‘  0  ) 

'  s  r  J 

d 

m 

0  -  0 
s  r 

0  '  0 

s  r 

3?ii  M  lA  a  a 
r  s  r  5 

2  2 
y  •  (o  -o  ) 

'  S  1 

2  2 
y  '  (0  '  0  ) 

'  s  r 

>  n  >  n 

r  I  S  S 
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The  trc.$sroiiclnf!on  as  ci  function  of  i  is  seen  to  bo  continuous,  in  TK.gnitude, 
evctywhptp  including  a  0.  At  ^  0,  botli  Equations  (B3)  arvi  f34)  irducp  to 


r  (0) 

rnr.ms 


4  2  2  2  2 

•32)t  M  M  a  a  y  '  (o  -  a  )  y  '  \o  ^  a  ) 

r  s  r  5  r  s  '  r  s 


Several  interestinq  special  cases  con  be  obtained  from  Equu’ion  (RS)  .  If  both 
mechonlcal  systenis  cte  undamped,  then  y  "  0  and 


I' 

*  m'  ,  ms 


t  t 
r  ! 


0  ,  a  ^  a 

I  s 


If  CT  a  ,  then 
r  s 


|-  (0) 

mt  ,  ms 


a  a  a 
r  5 


2tTy  F 


32n  M  M 


^  V  r  1 

Moo 

t  S  t  s  ' 


1  1 
1 - ?  I 

‘  4o  J 


T bus ,  i  f  o  o  and  v  0,  tlipn 
r  s 


I-  fO) 

mr  ,  rns 


- 

_  ^  o 


S  0 
r  5 


Equations  (86)  and  (88)  con  be  crjmbined  to  give  a  more  general  form,  namely 


'mr,ms^'^U<'  r'  ■  0 
r  s 


(0  ,  o  /  o 

6(o  -  o  )  \  ^  ^ 

^  '•  f 

’  00 ,  a  CT 

r  s 
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The  quantities  f  f  and  Cf  ore  ttie  bandwidths  of  the  two  systems  and  7  is  the 
r  r  s  s 

sum  of  the  bandwidths.  If  the  frequency  separation  jo^  ~  ^s|  niuch  larger  than 

7  ,  the  quntitl'y  a  above  will  be  much  smaller  than  tlic  quantity  b  above;  and 

the  quantity  c  nl.ove  will  ije  much  smaller  tfron  ttie  quantity  d  above.  Furtfier, 

2 

both  of  tlie  quantities  b  end  d  can  be  simplified  by  neglecting  y  so  that 


F' 


- - - - - 

16rt  M  M  CT  (o^  -  a  ) 
t  s  s  r  s 


F' 


m 


16n  a  ~  a  ^ ) 

ISIS  I 


In  this  cose,  Equations  (83)  and  (84)  are  appioxi  mo  tel  y 


(,  )  ^ 

ml  ,  ms 


A  O  ^  * 

16x1  M  M  (o  -  CT^) 
r  s  r  s 


-  F  '2ix  f  1 
s  s 
T  e 


sin  ?!T  a  T 


{90: 


2x,  F' 


1-  (-111) 
mr ,  ms 


4  ?  ? 

16x1  M  M  (a  -o') 
I  s  s  r 


T  P 


sIn2xxajT| 


wiieie  C  f  •  (  f 
r  r  s  s 


y  ^7 


u  -u 
r  $ 


I 


(bandwidths  do  not  oveilap) 


a  b 

c  d  . 


Note  that  the  largei  the  frequency  separation  jcr^  "  ''j  j  '  smaller  the  crosscorrelatlor 
(t)  (becomes,  cquotions  ( 90 )  ond  (91)  will  geneiolly  bx'  sufficiently  acc uro te 
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for  two  lightly  damped  systems  whose  bondwldths  do  not  overlop.  When  the  bond- 

widths  do  overlop,  that  is,  when  j  "  c',  |<  Y  *  dcmpings  of  the  Iwo  systems 

control  the  crossccrrelotion  V  (t)  of  thei.  responses  to  Impulse  excitations. 

mr  f  ms 

In  this  case,  the  more  general  Equations  (B3)  ond  (84)  must  be  used.  However,  when 
ja  “  ®  I  ^  )i  these  general  expressions  reduce  to  the  following: 


2n  F' 


I'  (t) 
*  mr ,  ms 


m 


32Tt  y  M  M  o  a 
r  s  r  s 


-  s  2r!f  T 

s  s  « 

e  cos  2n  0  T 

s 


(92) 


2rT  F 


r  (-1^1)^ 

mr ,  ms 


32n  y  M  M  o  o 
r  s  r  s 


cos  2it  0 


,hl 


(93) 


where , 

Icr  -  a  I y  (cj  '  a  )  ,  (bandwidtiis  overlap) 

I  !  s|  r  s 

b  «  a 
d  K  c 

t  C  f  •  C  f 
r  r  s  s 

Finally,  it  is  to  be  noted  that  the  slope  of  the  ctosscoirelntion 

continuous,  ingeneiol,  at  t  0.  Also,  since  for  i  ^  0,  the  frequency  of  the  function 

I  (t  )  is  cr  ,  and  -r  4O  the  frequency  of  the  func tion  I'  (-lTl)is  a  , 

mr,ms  s  mr,ms  r  ^ 

the  two  crosscorrelation  functions  may  be  quite  different  in  chorocter  .  These  characteristics 
of  the  ci  osscor  re  la  ti  on  function  ate  shown  in  Figure  8 
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F.  Crosscoi relation  of  the  Responses  of  Two  Single  Degree -of-Freedom  Systems 


to  Bandwidth  L'mitcd  V/hite  No'se 


Now,  considei  that  tlv?  same  two  independent  meclionical  single  degree  of  freedom 

systems  as  above  (Figure  7  ),  are  excited  by  bandwidth  limited  white  noise.  Then, 

cs  shown  in  Figure  9  ,  the  force  F  (t^,  input  to  the  first  system  results  in  an 

m 

acceleration  response  V*'  (f).  The  same  force  F  (t  *  i  ,),  input  to  the  second 

I  ml 

system,  with  time  delay  ^  ,,  results  in  on  acceleration  response  V/  (t) .  The 

I  s 

unit-impulse  response  functions  for  the  two  systems  ate  h  (t)  and  h  (t) 

mr  ms 


respecti  \'el  y 


If  the  first  lnf)ut  force  F  (t),  is  deterministic ,  then  the  second  input  force 

m 

F  (t  ‘  1  )  is  also  deterministic  and  has  tire  some  functionol  form  as  F  (t) . 

rn  I  m 

For  oxcifcplc,  if  F  (tt  A  sin  Prtff,  tlrcn  F  (t)  A  sin  2itf(t  *  t  Howu'ver, 
m  m  I 

if  F  (r)  reprt'senti  a  ranctom  signal  v.'hich  contains  a  continuous  band  of  frequencies 

over  a  finite  bondv.'id th ,  tlicn  the  composite  signal  F  '  t  )  over  all  frequencies 

m  1 

in  the  band  con  appear  Quite  different  from  F  (t).  Tc  obtain  F  (t  *  "r  ,)  from 

m  m  I 

F  (t),  each  ftequertcy  component  in  F  ft)  is  odvanced  in  time  by  a  constant  t  , 
m  m  l 

which  changes  tlie  p.iose  of  thof  frequency  component. 


The  amoijnt  of  phase  chonyo  2TTfT  ^  depends  upon  the  frequency  and  because 
tl-rere  ore  on  infinity  of  irequency  components  in  F  (t),  each  of  which  experiences 
G  different  phase  shift,  the  summation  of  oil  of  these  phase-shifted  freqrrency  com¬ 
ponents  will,  in  general,  produce  a  much  different  signal  than  the  original 
signal  For  small  values  of  t^,  the  phase  shi  ft  at  low  frequencies  wi  1 1 

be  small,  vdrereas  at  high  frequencies  the  shift  will  be  large.  Also,  the  larger  the 
time  Intel  vol  i  ^  the  iuigct  tlio  phase  change  for  all  frequencies  except  those 
iiovlny  periods  crppi  ox  I  mote  ly  equal  to  t  .  Tlnjs,  the  ctosscoi  relotl 


ion 
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1'  (t  *  T  )  between  F  (t)  and  F  (t  *  t  ,)  which  is  equivolent  to  the  auto- 
'  mm  1  m  m  , 

correlation  of  F  (t)  w! tli  time  delay  t  over  some  frequency  hand,  Is  greoter 
m  I 

when  the  band  contains  only  low  fiequencies  than  when  the  band  contains  only 
high  frequencies,  and  Is  greoter  for  small  values  of  i  ^  than  for  large  values  of 
T  ^  .  Also,  this  crosscorrelation  Is  dependent  upon  the  bandwidth  so  that  for  a 
given  i  and  a  given  center  frequency  of  tlie  band,  the  magnitude  of  the 
correlation  Is  greater  for  narrow  bonds  than  for  brood  botids  . 


The  similarities  of  the  tespor  ses  of  the  two  systems  In  Figure  9  reflect  the 

similarities  of  the  input.  Thus  the  crosscor relotion  I'  (t  )  of  the  responses 

ts 

V/  (l)  and  V/  (t)  is  greatest  when  tlie  crosscorrelation  I'  (t  ,  '  t  )  of  F  (t) 
t  s  mm  I  m 

and  F  ft  '  T  ,)  is  greatest.  Ho. /ever ,  ‘lie  crosscorrelation  P  ("r  )  is  also 
ml  rs 

governed  by  siml  la;  I  ties,  or  dissimilarities,  of  the  response  characteristics  of 

the  two  systems,  these  response  characteristics  being  the  resonant  freqirencies 

f  ,  f  and  the  dampings  i  ,  C  .  It  Is  expected  that  the  response  crosscorrelotion 
t  s  r  s 

f  (i  1  is  areotei  when  the  two  resononf  frequencies  f  ond  f  are  nearly 
'  rs  '  f  s 

equal  than  when  these  frequencies  ore  widely  seporated  .  The  effect  of  damping 
on  the  response  crosscorrelation  is  more  complex;  however,  Irr  general,  the 
response  correlation  decreases  os  the  damping  decreases,  since,  in  this  case,  the 
two  systems  become  mote  like  pure  tone  oscillators  with  zero  bondwidth,  for 
which  tlie  crosscorrelotion  is  zero  If  the  tv.'O  resonant  frequencies  are  different. 


Equation  (76),  tire  general  expression  for  the  system  shown  in  Figure  6  ,  can 

be  applied  to  the  syste  m  being  considered  here  by  using  only  the  first  term  which 

is  the  crosscorrelation  of  two  single  degree  of  freedom  systems  to  the  same  Input 

force  F  (l),  hence 
m 


I  (0  (2rTf) 


/  I  l-r  I  ( 

y  mf,ms  /  m 


T 


T  2^ 


2 


(94) 
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The  crosscorrelation  ['  (t  *)  of  the  unit  impulse  response  functiors  is  obtained 

mr  f  ms  i 

from  Equations  (83)  and  (34)  by  setting  F  “1, 

rn 

-  F  2rrf  T  2 

‘mr,ms^"2'  "  ^  *  [a  cos  2n  t  ^  ^  b  sin  2  } "  2 ^  ° 

.  -f  2iTf^T  111 

'mr,ms^‘  I"'  2r  ®  '  '  |c  cos  2n«T  Jt  2I  ^  d  sin  2n<Tj  t  2]  J  ,  t  2^  0  (99) 

where  a,  b,  c,  and  d  ore  os  defined  for  Equations  (33)  ond  (34)  with  F  "  1  . 


Substituting  Equations  (97),  (98)  and  (99)  Into  Equotlons  '96)  gives 


r  (3)  -  (2nf) 
rs 


00 

■'m/ 


-  f  2itf  7  $m  2itf  (t  -7  -)  -  Sin  2Kf,  (7,  -  7  _) , 

ss4„  u12  L12d7_ 

e  cos  /nff  7  — -s -  2 

‘  ^  1  -  ’2’ 


/CD  -  L?»rf  r  sin  2nf  (7  ,  -  7  _)  .  -  ,  ,  . 

'  o  1  2-sin2nf,  (7,-7^), 

=  i.n  ^ - _jJ_J - iU-r  j 

<’  1  •  ’  2' 


00  -  f  27T  f  7  - 

'  f  e  '  'I  2 

m  J  cos  /TTfT  ^  M  2 

o 


sin  2nf  (7  ,  -  7  ^)  -  sin  ?nf  (7  ,  -7 
_ u  I  2 _ L  I  2 

<’  1  ■  ’  2’ 


d7 


>  d  5  /  e  ■ 2.,  1  ,  I  I  -  -  2>  -  2>  d,  ^ 

^  Jr.  '  2  ,  .  ^ 

o  •  '  (7,-7^) 


1  2’ 


(100) 
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Spectral  Density  of  Response  Differences 


A  convenient  means  for  presenting  flie  frequency  response  transfer  relationships 

existing  between  two  acceleration  response's  W^(t)  and  V/^(t)  at  poltits  y^  and  y^ 

on  a  complex  structure  is  fc  und  by  computing  the  spectral  density  S  ,  v(f)  of 

«  •  •  •  *  ' 

the  difference  between  V/  (t)  and  W  (t).  Equations  (42)  and  (43)  define  the 

autocorrelation  5'  y)  and  the  spectral  density  S,  .(f)  of  the  response 
(r  -  s)  '  (t  -  s) 

•  •  •  • 

difference  V/  (t)  -  V-'  (f).  Expond  Equc.tion  (42)  and  obtain  the  arjtocorrelotion 
(  s 


1,  ,fi  )  [wff)-W(t)  fw  (t  r  r)  -  W  (f  4  r)| 
(i  -  s)  r  s  J  1  r  s  J 


-71 - n 


w  (f)  w  (t  ■  7  )  -  v;  (t)  v;  (r  .  T  )  -  w  (t)w  (t  <  i  ) 

(  t  r  s  SI 


W  (t)  V,'  (f  •  1  ) 

s  s 


f.  (t  )  r  (t  )  -  r  (i  )  -!■  (t  )  '  I-  (-.  ) 
(r  -  5  r  rs  sr  S 


(101) 


The  spectral  riensity  ^  ^^ff)  of  ‘he  response  difference,  by  inspection,  is  written 


as 


-(■  - 


CD 


/ '  - 


-i  2ii  fi 
)  e  di 


cn 


r  IS 


sr 


(102) 
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wlosc  exomination  of  EquoMon  (102)  revcaU  a  useful  tool  for  presenting  meosured 
esponse  data . 


For  example,  consider  tv/o  independent,  but  identical  systems.  V/hen  two 

independent  random  input  forces  F^(t)  and  t’^ft)  are  applied  to  the  systems,  their 

•  •  •  • 

acceleration  respon'^es  W  (t)  and  VJ  (t)  are  uncor related  .  For  this  case,  the 

r  s 

••  •• 

spectral  density  S  ^  ^^(f)  of  the  lesfsonse  difference  'A'  (t)  -  W  (t)  becomes 

simply  the  sum  of  the  separate  spectral  densities 


S,  .(f)  S  ff)  '  S  (f),  (unconelated), 
(r  -  s)  r  s 


(102) 


which  is  what  would  be  expected,  since  uncoircloted  responses  ore  Independent 
of  their  relative  phose,  and  hence,  tl^e  sign  of  theli  combination;  that  is 


S,  .(f)  S  (f)  *  S  (f),  (uncorr.elored) . 
(r  t.s)  t  «  ' 


For  the  case  of  perfectly  correlated  input  forces  v/Ith  identical  Input  spectral 

densities  S.(f)  '  S_(f)  S(f),  the  cross-spectral  densities  S  (f)  and  5  (f) 

I  ^  rs  sr 

and  the  spectral  densities  5  (f)  and  S  (f)  ate  eqi.-al,  hence  tire  spectra!  density 

r  s  ' 

5^  jff)  of  the  response  difference  is  zero. 


s.  .(f)  S  ff)  -  S  (f)  -  S  (f)  -  S  (f) 

(i  -  s)  r  rs  sr  $ 


S,  .  ■  0,  (correlated) 
(r  -  s) 


(104) 
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The  response  difference  spectral  density  5^^  is  sensitive  to  the  relative 

correlation  existing  between  the  frequency  spectra  of  the  two  resfxrnses  W^(t) 

and  W  (t).  Wiien  5,  .(f)  is  large  at  some  ftcq  jency  f,  there  is  o  lack  of 

s  (i  -  s) 

correlation  or,  rather , "transmission  loss"  of  that  particular  frequency  com¬ 
ponent  between  the  points  y*  an  I  .  Tire  distribution  of  various  frequency 
components  of  vibration  can  be  traced  fhroirghout  a  conrplex  structure  via  this 
method,  which  would  result  in  tiie  location  and  identification  of  excita’ion  sources 
from  data  measured  at  selected  points  on  a  structure.  Figure  lO  schematically 
illustrates  tire  utility  of  this  method. 


The  spectral  density  S,  ,(fi  of  the  respemse  difference  W  (♦)  -  W  (t)  is 

indicative  uf  the  'transmission  loss"  of  the  stiijcPital  segmerri  between  the  response 

points  y*  and  .  Dividing  Eejuation  (10?)  by  Eqtjntion  (103),  whicli  is  the 

maximuei  value  of  S,  ,ff),  yields  ci  nondimensional  coefficient  TL,  . 

(r  -  s)  '  (r  -  s) 


TL 


fr 


-s) 


^1l: 

S  ff) 


.(f) 


S  (f) 
s 


(105) 


which  evoluotes 


the  transmission 


and 


y  • 
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DISTRIBUTED  PRESSURE  EXCITATION  AT  DISCRETE  FREOUENCiES 


Consider  now  the  case  of  a  d'stributed  pressure  looding  which  has  a  discrete  frequency 
spectrum  end  let  P|^(x,t)  denote  the  component  of  the  pressure  opplied  at  the  point 
"x  at  frequency  fj^  .  An  expression  for  the  acceleration  response  W^(t)  at  a  point 
can  be  obtained  from  equation  (4)  by  replacing  the  point  force  by  its 

equivalent  P[^(x/  0  '  bA  (x)  ,  where  dA  (x  )  denotes  an  incremental  area  over 

which  the  pressure  acts,  and  by  replacing  the  summation  over  n  by  an  Integration 


over  area,  giving 


vy  (t) 

r 


f,  )  ■  Pi  (x  ,  t)  •  dA  (  X  )  (106) 


In  Equation  (106),  the  areo  A  denotes  tfie  entire  area  of  tire  veii'cle  over  whiclr  the 
fluctijotirrg  pressures  act.  For  purposes  of  developing  err.pirlcal  correlation  techniques, 
it  is  convenient  to  break  up  this  surface  area  Into  zones  or  areo  segments  A^  as  shown 
in  Figure  II  .  The  choice  of  the  size  of  these  zones  will  generally  depend  upon  the 
surface  pressure  correlation  lengths,  crosscorrelations  of  tlie  responses  along  different 
structured  transmission  paths,  exciratlon  frequ-ency  ond  bandwidth.  Discussions  on  these 
points  ate  presented  In  Section  IV.  Tire  Integrol  over  the  vehicle  areo  A  can  be 
written  as  a  summation  of  Integrols  over  each  of  the  area  segments  A^,  and  thus  the 
expression  for  the  accelerotion  response  V/  (t)  iii  Equotion  (106)  can  be  written  in 
the  form 


;  f^)  P^('1  ,  •)  •  dA  ("^  )  (107) 


The  pressure  ,  t)  can  be  written  in  complex  form  in  a  me  rner  similar  to  that  used 

for  the  cu  iplex  force  F  ,  (t)  given  in  Equation  (19)  , 
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(^,  t) 

Also,  tlie  complex  fotm  of  o  ( "?  , 
given  by  Equation  (7)  ,  namely 


j  i*  (  X  )  ^ 

)|  •  e  e 

con  be  written  in  a  form  similar  to  that 


(1^ 


'i.' 


i9( 


y.' 


(109) 


Substituting  Equation  (103)  <  nd  (10.^)  into  Equolion  (107)  gives  tlie  following 
expression  for  the  acfeleta*ion  W  (t), 


W  (t) 


Equation  (110)  is  seen  to  hove  tbe  same  form  as  Eguotion  (20)  and  tlujs  t!ie  expression 
for  the  meon-squore  occeletation  response  V/^^(t)  is  similar  to  that  given  by 
Equation  (22).  Thus  for  the  case  of  a  distributed  pressure  looding, 


W^(f) 


+EE// 


o  (?,■?■ 


A  A 

n  n 


EEE  // 

k  n  m  / 

(n  /  m)  n  rn 


'  ,  y^  ;  f|^)  dA(  x  )  o'A  (  x'  ) 

(111) 

O  d,.  (■?  )  dA(l^  ) 
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where 


0(x  ."y^  ;  1^)  |a  1?  ;  f|^)  |  •  | a  (5?'  ,  ^  j  C^)  j  •  |  P^(^'  ) 


•  cos 


6("^  '  '■  ~  '"^r  '  ^  ^  '  'f’k  ^ 


Similarly,  the  crosscorrelation  I'  (r)  of  the  accelerotion  responses  W  (t),  W  (t) 

rs  r  s 

at  y^,  respectively  lies  the  some  form  os  Equation  (29)  and  is 


-\\  (^)|  |Pi,('^'  )i  • 

•  cos  f^)  -  0("J',  7^;  -  4,^  (oj  (11^) 

It  is  to  be  noted  that  Equation  (113)  ond  (114)  reduce  to  Equation  (111)  and  (112) 

when  y^  -  y^.  Cor, siderlfig  the  ger.ernlity  of  Equations  (1 10)  -  (1 14)  for  the 

mean-square  response  and  crosscorrelation  of  response,  these  expressions  ore 
quite  compacr;  and  being  expressed  in  real  form,  they  ore  easily  interpreted. 

Normal  Mode  Rc-prcsontatlon  of  Spectral  Density  and  Cio.sF^correlatlon 

When  equations  of  this  type  are  presented  in  the  literotuie,  they  are  usually  presented 
in  terms  of  normal  modes  of  the  structure.  To  show  hov/  these  equations  can  be  written 


in  terms  of  normal  modos,  it  is  convenient  to  v/rite  the  expression  (114)  for 

Q(  X  ,  x’  ,  y  ,  /  an  aherna?c  for,-n.  In  order  to  shorten  the  notation, 

f  s  k  ^ 

let  a  ,  a*  ,0  ,0*  he  used  to  denote  a  (x  ,  y  ;  f,  )  ,  a  (  x*  ,  y  ;  f,  )  0(x  ;  y  ;  h  ), 

r^srs  rk  $k'  rk* 

and  0(x'  ,  y^  ;  f^ )  respectively.  Nd  te  first  that  the  real  and  Imoglnary  forms  of 

a  and  o'  ore 
r  s 


^l",}  Kh'”®, 

Jkl  kl'i"8; 


(115) 


Expond. ng  tl.e  cosine  function  in  Equc'ion  (114)  ond  using  Equation  (115)  to  simplify 
the  resulting  expression,  gives 


or”'  i  \ 

Ql>  .  -  .  y,  .  y,  .1^1 


.“.I'^Kl]  '  j '’i,* '  •  I  ■  cos  »^( ») )j 


I  Sin  <?>|^  (  X  )  -  (x 


As  shown  in  n  previous  exanyple  for  a  uniform  beom,  the  form  of  the  tronsfer  fu  iction, 
when  expxTnded  in  teirr..  -.f  ommol  modes  of  the  structure,  Is  that  given  by  Equjtlon  (13) 
Thus,  foi  a  genrrni  stiuctum,  r,  can  be  written  as  follows 


a  "  ( ?rr  f  ) 
r  k 


,?  "A’ 


/  .  ~r~ 


(117) 
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and  the  comparable  equation  for  Is 


a  -(2ttfj^)  2L  “T 
M  l  M 


(113) 


The  subscripts  *'  and  p  denote  the  mode  number,  and  k  tlie  generalized 

modal  stiffnesses,  t\i  ^  C")  and  0  (  X  )  the  mode  shapes.  The  complex  fuiiction 

H 

H„  (f)  is  defined  by  Equotion  MOl  /  and  the  real  and  imaginary  parts  of  H,,  (fj^)  ore 

'-ft)  K^'T) 


(119) 


where  the  function  j  Hj,  (flj  is  defined  by  Equation  fll),  Similotly  for  the 
function  (f|^) .  Substi  tutlng  Equations  (117)  -  (119)  into  Equation  (116)  gives 
exfxinsion  of  0(  x  ,  x  '  ’  ^k^  terms  of  normol  modes  of  the  structure 


OP  oo 


-4  -> 

0(  X  ,  X-  ,7  ,7  ;  f,  1 

t  s  k 


<^'k' 


;•  I  n 


y{',l  X 

ip  ,p 


p,(7)Mp,(t')| 


(120) 


where  the  functions  X  and  *X  are  defined  os  follows: 

‘e  1'^ 


tV„  ■  K'ffc'f  •  l^('k>r  ■ 


(121) 
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(x  )  -<1> 


-4, 


) 


(r?2) 


It  !s  to  be  noted  that  when  ^  ,  the  second  term  in  Equation  122  is  equal 

to  zero .  By  using  Equations  (11?)  -  (1  22), Equations  (111)  and  (113),  the  mean- 
squaie  response  at  a  point  on  the  stnjcture  and  the  ciosscotielatlon  of  responses 
between  tv.'o  points  on  the  structure , can  f)e  expressed  In  terms  of  normal  modes 
of  the  stnrctuie  . 


DISTRIBUTED  RANDOM  EXCITATION 

The  analysis  de.eh)p*--d  ohove  for  distributed  pressures  ol  discrete  frequencies  can 

be  exterrded  to  the  cose  of  distributed  randomly  varying  pressures  acting  over  tlie 

whole  surface  area  A  of  a  complex  structure.  The  precedlfig  section  on  point 

force  random  excitation  contains  all  of  tlie  necessary  development  for  determinrng 

the  important  qrrantilles  and  parameters  necessary  for  describing  the  response  of 

a  structure  to  distrib'rted  loadings.  To  accornpllsli  this  ti  insltlon,  replace  the 

poin*  random  forces  f  ft)  which  ate  applied  at  the  prjints  x'^  by  therr  equrvalent 

pressures  P(x*  t )  d  A  ("j^) ,  wliere  d  A  (x^  denotes  an  i  ncremen  to  I  at  ea  over  which  the 

prP5<,ij,f.  acts.  The  surface  area  A  is  represented  by  xoncs  or  area  segment  A^^ . 

Thus,  the  equation  for  tlie  accelerotion  response  V<'^(t)  at  o  point  y^  can  be 

obtoined  from  a  form  of  Equation  (34)  by  on  integration  over  each  area  segment 

A  ond,  tlren,  a  summation  over  the  N  area  segments, 
n 


W^(t)  X  f  h(x‘,y*;a)  Pfx,t-'T)dT  dA(x).  (123) 

ri  I  A  •'O 
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The  autocor'olation  function  I'(t  )  of  the  response  W^(t)  con  he  defined  In  a 
manner  siniilnr  to  the  dovelopmr  nt  of  Kfjuntion  (5S),  hence,  the  equation 


[•(■r  )  =  W^(t)  W^(t  4  T  ) 


N  N 

m  ^  1  n  1 


2 


) 


•  P(x ,  t)  P(x  ' ,  t  4  T  ^  -  1  2  *  T  )  di  ^dA(x)dA(x  ) .  (124) 


And,  the  spectral  density 


S  (f)  of  W  (t)  Is  defined, 
r  r 


from  Equation  (33),  os  the 


equation 


S  (f)  -  2 


r  (T  )  e 


-i2nfT 

di 


N  N  r  r 

E  r// 

-  1  n  -  1  -A  -A 

n  n 


(2nf) 


i  2nfT 


r  )e 


dr 


-  i2rTfT 


h(x 


dT 


OO  ^  -I2ltf(T  -  ^  '  T  ) 

.2/  P(x,t)  pfx',  t  *  T  ,  -  T  2  *  T  )  e  d(T  ^-7  2  +  T  ) 

-  CD 

I 

.  dA(x )  dA^  ) 


6/ 


N  N 

t 

m  _  1  n  1  A  A 
n  m 


If 


a  (x,  y  ;f)  a  (x  ,  y  ;f)  •  S(x,x  ,f)  dA(x)  dA(x  ) 


Thus,  flip  nieon-squatc  occelerotlon  rpspo^^e  becomes 

- r: -  CD 


W  (f)  I  S 

'  {  ' 


(f)df 


026) 


Equations  (1  24)  and  (1  25)  aie  general  expressions  for  tlie  autocorrelation  l^(  t  ) 
ond  spectral  density  S  (f)  of  the  response  W^(t)  of  o  point  y^  due  to  N  distributed 
loadings  over  the  area  segments  A^  • 

Equation  (125)  can  be  written  in  terms  of  the  normal  mode  shapes  (})  (x)  of  the 

— ♦  ■  -^1 

entire  complex  structure.  In  this  cose  the  tronsfer  function  a  (x,  y^;  f)  can  be 
written  os 


a  (x,  y  ;  f) 


00 

E 

q  -  I 


or  (f)  (y  ) 't'  (x), 
qr  ^q  't  q 


where  a  (f)  is  the  generalized  transfer  function  of  the  q  -  th  mode,  (|)  (y  ) 

is  the  q-th  modal  deflection  of  y*  and  (t  (x)  is  the  q-th  normal  mode  of  the 

r  q 

structure.  Hence,  Equation  (125)  becomes 


00  00  N  N 


S  (f) 


E  E  E  E 

q  1  j  "  1  m  1  n  1 


.  (f)  . 


,  ^  r 

<t>  )  <1>.  (y  )  /*  /*  S(x,  X  ;f)  (|)  (x)  ipjx) 

^  \  ^  J  J  q  I 


A  A 
n  m 


*  dA(x )  dA(x  )  . 


(127) 


(125) 
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EquoHon  (129)  can  also  be  written  In  terms  of  the  normal  mode  shapes 
Ip  (x)  as  the  equation 

q 

CD  oo  N  N 

-  E  E  E  E”V" 

q  "  1  j  “  1  m  ”  1  n  -  1 


(|»q(y^)  (fij  (yj  J  J  S{*x,  ^  ;  f)  <})  ^  (x)  (}).{x  )  .  dA(x)  dA(x  )  (130) 


A  A 
n  m 


Response  of  a  Beom  to  Disttibut^d  Loading 

Consider  a  uniform  beam  of  length  L,  excited  by  u  Jistributed  pressure  loading 
of  Intensity  P(x,f)  per  unit  length.  The  acceleration  response  W^(t)  at  a  point 
X  on  the  beam,  wMcIi  Is  in  its  first  normal  mode  (f^{x)  only,  subject  to  the 
distributed  looding  P(x,t)dL(x),  is  written 


W^(t) 


V  „..2  '7- 


i  t-r)dr  dL(x) 


(131) 


The  spectral  density  S,(f)  of  the  response  in  the  first  mode  from  Fquation  (127), 


S^(f)  '  ja^^ff)]  ^  <p|  (^^)  J  ^  S(x,  X  ;f)  (})^  (x)  (x  )  dL(x)  dL(x  '  032) 


where  Sfx^  x  ;f)  is  the  cross-spectrol  density  of  the  pressures  at  T  and  x  , 
a^ff)  is  the  generalized  transfer  function  of  the  first  mode,  (x^)  is  the  modal 
deflection  at  x*  and  (}>j(x*)  is  the  first  normal  mode  sliape  of  the  beam. 
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Now,  the  double  integral  over  the  length  L  of  the  beam  from  Reference  2 
is  of  the  form 


Sp(f)L^i,^^f)  -  J  (  s(ir,t ';f)^^^),5>^r;')  di^dK^')  (133) 

0  *0 

2 

where  S^(f)  is  the  spectrol  density  of  the  pressures  and  (f)  is  the  'joint 
acceptance  squared  '  of  the  first  beam  mode  ^^(x)  subjected  to  the  distributed 
pressures  P(x,t)  .  Therefore,  substitution  of  Equation  (133)  into  Equation  (132) 


L  L 


gives 


S/0  -  <^/(",)  |oj(f)|^  S  (f)  j^,^(r)  • 


(134) 
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SECTION  IV 


GENERAL  METHODS  FOR  EM  PI  RICAL  CO  RI^E  lA  T I O  NS_ 

^ACKX.ROUNp 

Ttie  piocf'5'i  of  empirical  coiielotion  of  viLratlon  or  other  data  generally  begins  with 
the  comparison  of  the  data  found  for  specified  conditions  on  one  or  more  flight  vehicles. 
If  these  data  show  any  degree  of  similoiity,  or  can  be  normalized  to  show  similarities, 
the  engineer  is  encouraged  to  utilize  them  as  o  predictive  tool  for  a  future  vehicle. 
Correlations  of  this  type  have  been  made  throrjghout  engineering  history,  and  though 
a  rational  basis  for  their  use  is  not  always  easily  found,  they  have  often  been  extremely 
useful  and  surprisingly  accurate  .  As  a  result,  one  might  question  the  real  reasons  behind 
the  oppi  i  cohi  I  i  ty  of  the  empirical  correlation  technique  to  mony  dlvt.’rse  vehicles. 

In  the  exomlnation  of  doto  from  sirrgle  engine  figlrter  aircraft,  multi-engine  bombers, 
and  launch  vehicles  for  space  flight,  Reference  3  ,  one  con  find  many  similarities 

despite  the  gross  apfrorerrt  differences  In  the  structural  types.  The  existence  of  tirese 
similarities  suggests  that  many  slnrcturnl  types,  even  thorjgh  rodlcoily  different  in 
geometry  and  other  opporent  design  foctors,  do  possess  certain  properties  wtrlch  Irave 
statistical  similarity.  Petliaps  the  most  readily  identifiable  of  these  properties  is  the 
high-frequency  Input  impedance  of  complex  structures  which  tends  to  be  purely  resistive 
as  the  structrrre  begins  to  behove  ns  n  transmission  line.  Preliminary  investigations 
srrggesf  tlrat  the  application  of  the  model  scale  laws,  modified  by  consideration  of 
the  volumetric  density,  will  enable  direct  compotison  of  doto  over  a  wide  range  of 
velricle  types.  In  fact,  cempor'son  of  the  input  impedonces  of  a  submarine  and  Snark 
ft/selo.]e  oppeor  to  correlate  well  wheti  adjrjsted  for  these  deterministic  parameters. 

It  is  impo'tant  to  recogr.ize  that  the  validity  of  an  empirical  correlation,  porticularly 
when  the  correlo  Ion  Is  to  be  utilized  as  o  predictive  tool  for  a  future  vehicle,  is 
entirely  dependent  on  the  statistical  similarity  of  those  structural  properties  which  are 
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significant  to  vehicle  vibration.  Therefore,  any  correlation  of  structnrol  vibration 
which  accounts  for  the  gross  structural  porometers,  at  least  in  terms  of  the  model 
scale  laws,  volumetric  density  and  dynamic  magnification  factor,  is  applicable  to 
other  vehicles  as  long  as  their  structures  actually  exhibit  similar  properties.  Thus, 
empirical  correlations  may  often  be  applied  to  g  greet  variety  of  geometrical 
configurations . 

The  next  subsections  will  review  two  of  the  simpler  correlations  which  have  Ix'cn 
made  in  the  past  and  are  follo.ved  by  the  development  of  a  generol  method  for 
empiricol  correlation. 

'ie 

Direct  Correlation  of  Noise  and  Adjacent  Structural  Vibration  in  Aircroft. 

One  opprooch  to  the  correlation  of  vibration  response  of  aircraft  v^Ith  external 
noise  has  been  suggested  previously  by  Convair  In  References  4  -  6  .  Tlris 
correlation  was  obtained  by  direct  comporison  of  the  exlernol  sound  pressure  level 
on  the  B-58  resulting  from  jet  engine  noise  with  the  vibratory  response  of  adjacent 
internal  structure.  Tlius,  vibration  meosurements  in  the  nose  of  tlie  oircraft  were 
compared  directly  to  sound  pressure  levels  measured  on  the  external  skin  of  the  nose 
section,  and  vibration  measurements  in  the  aft  end  of  the  aircraft  were  compared  to 
sound  pressure  levels  measured  on  tire  external  skin  of  the  oft  section,  etc.  Since 
the  external  noise  environment  on  the  B-53  and  other  types  of  jet  aircraft  varies 
considerobly  from  the  relatively  low  noise  levels  forward  to  the  very  high  levels 
toward  tfie  aft  end,  a  comparisjn  of  this  type  includes  a  range  of  external  noise 
levels  of  the  order  of  approximately  30  db . 

A  similar  comparison  of  external  noise  and  internal  vibration  from  Snark  data  in 
References  7  and  8  Is  given  for  two  frequency  bands  in  Figure  12  .  The  lower 

pressure  levels  ernd  their  assocloted  vibration  data  are  obtained  at  forward  fuselage 

*  Taken  from  i’.efprence  3. 
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stations,  and  the  higher  sound  pressure  levels  and  the  ossocioted  vibration  data  are 
obtained  at  oft  fuselage  stations  adjacent  to  the  rocket  booster  exhaust.  It  should 
be  noted  tiiat  the  small  variations  in  sound  pressure  level  at  one  microphone  p3sItion 
result  from  data  of  repeated  firings  ortd  that  tlie  acceleration  data  Include  some  from 
transducers  attached  to  airplane  structure  and  oriented  In  various  directions. 

As  might  be  expected,  the  dato  exhibit  consideroble  scatter,  so  that  it  Is  possible 
only  to  estimate  a  rend  n-  through  the  points;  however,  more  significant  tiran  the 
scatter,  is  the  slof  trend  line.  Note  that  if  the  acceleration  amplitude 

increases  by  a  fa  or  c  for  on  increase  of  20  db  spl ,  a  direct  linear  relationship 
exists  between  the  internul  strucfurol  vibrotion  amplitude  and  the  adjacent  external 
sound  pressure  amplitude.  In  this  event,  the  trend  line  through  the  data  gives  a 
constant  of  proportionality  between  external  sound  and  internal  vibration.  However, 

If  the  trend  line  has  a  slope  of  less  than  unity,  as  Is  the  case  of  the  low-frequency 
data  of  Figure  12  ,  a  direct  linear  relationship  between  internal  vibration  and  adjacent 
externol  sound  pressure  cannot  be  proved  by  the  correlation. 

Severoi  factors  might,  individuolly  or  collectively,  he  responsible  for  slopes  less 
than  unity.  These  factors  include: 

1  .  Structure  borne  transmission  of  vibration  from  high  external  noise  level  oreos 
to  those  of  lower  external  noise  level; 

2.  Nonuniform  structure  throughout  the  fuselage  (which  undoubtedly  contributes 
to  scatter ); 

3.  Non  linear  response  of  tire  structure; 

4.  Unfortunate  selection  of  t.ansducer  location  with  respect  to  model  response. 

Before  considering  any  of  these  factors  in  more  detail,  it  Is  helpful  to  examine  the 
slopes  of  the  trend  curves  as  a  function  of  frequency.  Figure  13  summarizes  the 
results  of  this  type  of  direct  correlation  for  four  aircraft  type  structures,  Includinq 
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the  B-66  of  Reference  9  ,  together  with  Snork,  B-58  and  B-52.  The  upper  portion 

of  the  figure  gives  the  slopes  of  tire  trend  curves  and  the  lower  portion  gives  the 
average  rms  occeleratlon  for  each  octave  band  when  the  adjacent  external  octave 
band  spl  Is  150db.  It  is  clear  from  Figure  13  that  some  correlation  betv^een  adjacent 
external  noise  may  exist  above  150  cps  for  the  Snark  and  at  the  higher  frequencies  for 
the  other  aircraft  where  the  slopes  of  the  trend  curves  approach  unity.  Since  a  large 
number  of  randomly  selected  transducer  locations  ore  included  In  the  various  surveys, 
it  is  doubtful  that  the  occidentol  location  of  transducers  with  respect  to  the  various 
vibrot'on  modes  accounts  for  the  Iov;er  slopes  at  the  lower  frequencies.  Similarly, 
since  several  experiments  give  the  some  general  conclusion  with  regard  to  the  trend 
curves,  it  is  not  felt  that  structural  nonunI form! ty  is  a  particular  foctor  in  the  deter¬ 
mination  of  the  slope  of  the  trend  curve,  altirouyh  it  undoubtedly  Is  a  most  significant 
factor  In  the  scatter  of  tlie  Individual  data. 

The  role  of  nonlinear  behavior  In  the  structure  cannot  be  evaluoted  with  respect 
to  the  trend  curves  from  the  available  data;  however,  it  is  considered  to  be  less 
important  than  the  fact  that  the  response  ot  any  generol  position  in  the  structure  Is 
given  by  the  sum  of  the  noise  energy  transmlf'ed  directly  to  adjacent  structure,  plus 
the  noise  energy  received  at  more  remote  locations  and  transmitted  as  vibrotory 
energy  through  the  structure  to  the  position  as  illustrated  irt  Equa tion  ( 1  dT,)  .  Thus, 
at  higher  freqi;encles  better  correlation  would  be  expected  between  adjocent  ocoustical 
excitation  and  respijnse  liecause  the  vibrational  energy  transmitted  through  the  srructure 
from  remote  locations  has  been  signifierrriti  y  attenuated  and  mokes  cmly  a  minor  con¬ 
tribution  Conver.sely ,  ert  low  frequencies  tfie  vibrational  energy  transmitted  frorn  the 
areas  which  have  the  highest  external  noise  levels  would  be  expected  to  exceed  the 
local  excitation.  This  generol  result  might  be  anticipated  since  the  attenuation  of 
vibratory  energy  transmitted  by  bending  waves  otong  the  fuseloge  is  essentially  con¬ 
stant  per  wave-length.  Thus,  tire  high-frequency  energy  suffers  considerably  riiOre 
attenuation  than  does  low-frequency  energy  v/hen  both  are  transmitted  for  the  sam.e 
distance  throirqh  tlie  fuselage. 
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Coitelafion  of  Ballistic  Missile  Launch  NoUe  end  Response 

The  vibrotion  dafo  given  in  Reference  3  for  structuie  located  in  the  forward  region  of 
bollistic  missiles  during  botli  lounch  and  maximum  dynamic  pressure  flight,  consistently 
showed  that  the  vibration  environment  of  the  large  end  heavy  missiles  was  significantly 
less  than  the  vibration  of  the  smaUer  and  lighter  missiles.  This  suggests  that  the  next 
generotlon  of  space  vehicle  launch  platforms,  which  are  anticipated  to  be  considerably 
larger  tfion  tire  picsent  ICBMs,  might  have  even  less  severe  vibration.  On  the  other 
bond,  chonges  i'l  lounch  configuration,  engine  parameters,  or  maximum  q  could 
result  In  Increases.  It  Is,  therefore,  desiroble  to  determine  on  empirical  relationship 
iiptween  tfie  response  in  fotv/ard  structure  and  the  forcing  function,  appropriately 
modified  by  r''e  vehicle's  structural  parometers. 

It  is  believed  that  an  empirical  correlation  can  be  found  by  equoting  the  available 
power  of  the  forcing  function  wittr  the  mechanical  energy  stored  In  the  vibratory  system 
and  pov^cr  di.'-'.sii'rdnd  tfrrough  domping.  It  is  fundomontal  that  these  quontities  ore 
related  since  all  of  the  incident  acousticol  power  which  is  absorbed  by  the  vehicle 
must  result  In  vibration  and  be  dissipated  through  damping. 

r  r  tills  special  cose,  an  alternative  opprooch  to  the  correlation  of  the  launch  acoustical 
forcing  function  wlt!i  vehicle  response  oppears  to  give  encouraging  results.  It  is  well 
known  tbot  corrplicated  (and  even  simple)  structure  has  many  resonances,  distributed 
throughout  the  frequency  range,  which  occur  at  frequencies  above  the  fundamental 
resonant  frequency.  Consequently,  tlie  response  of  a  structure  to  a  random  forcing 
function,  which  has  energy  approximately  equally  distributed  over  a  brood  frequency 
ronge,  exhibits  man^  resonont  peaks.  Similarly,  the  majority  of  the  missile  launch 
data  In  Figure  14  are  representative  of  the  vibrotion  amplitudes  associated  with  one 
or  more  of  the  many  resonances  which  occur  in  the  missile  . 

''r  Taken  fror;  Keferenc’  1. 
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It  is  well  known  that  the  natural  vibration  characteristics  of  many  complex  structures 
can  be  approximoted  by  individual  consideration  of  each  resonance  or  mode  of 
vibrotlon,  ossuming  It  to  be  essentially  unaffected  by,  or  decoupled  from,  ony  other 
mode.  The  response  of  the  vcl'.Icle  in  any  one  of  these  modes  con  be  obtained  from 
expressions  similar  to  tliose  developed  for  the  single  degree  of  freedom  if  an  appropriate 
definition  con  be  obtained  for  the  amount  of  the  total  moss  which  Is  Involved  In  actuol 
vibratory  motion,  and  for  the  effective  or  "generalized  force"  on  the  vehicle.  Viewed 
in  this  perspective,  th»>  totol  response  of  the  vehicle  ot  any  location  is  simply  the  sum 
of  the  contributions  from  all  of  the  vehicle's  vibrotory  modes. 

The  mean  square  displacement  response  of  a  single  degree  of  freedom  system  to  a 
continuous  random  foicing  function  Is  given  by; 

“I  ^“n  7(f? 

"  ^  - 

where 


and 

The  mean  square  acceleiallon  ratio  (a/g)^  associated  with  this  response  can  be  easily 
obtained  from  the  Erjnation  .TJrove: 


2 

X  is  the  mean  square  displcicement 
Q  is  one  divided  by  twice  the  domping  ratio 
u  is  the  natural  frequency  in  radians/second 


F(f)  is  the  mean  square  force  per  cycle/second 


k  is  the  stiffness  of  the  system. 
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whei  e 


W  is  the  weight  of  the  moss. 

In  order  to  apply  this  simple  concept  of  the  single  degree  of  freedom  system  to  the 
empirical  correlation  of  multimode  missile  vibration,  it  Is  necessary  to  consider  the 
relationship  between  )he  generalized  force  and  generalized  mass  of  the  vehicle  with 
the  single  forcing  function  and  single  mass  in  the  single  degree  of  freedom  solution. 

For  example,  the  genetollzed  force  on  a  panel  in  o  space  vehicle  resulting  from  acoustic 
excitation  depends.  In  addition  to  the  actual  magnitude  of  the  sound  pressure,  primarily 
upon  the  spatial  corioloflon  of  tlie  sound  pressure,  c  tlie  distances  over  which  the 
pressure  is  in  phase.  In  relationship  to  the  distorice  betv/een  modes  In  fhie  \  el's  bending 
response  .  It  Is  also  hnown  that  when  acoustic  phcnomeno  are  similar,  except  '‘or  a 
scale  factor,  their  spatial  correlations  ore  similar  when  compered  on  a  wave  number 
basis.  Here,  the  wave  number  (kr)  is  equal  to  the  number  of  rodions  per  unit  distance 
(rV c  "  k)  times  a  distance  which  is  chosen  to  be  clsaractcristic  of  the  size  of  the  object 
or  missile.  For  the  purpose  of  this  correlation,  r  is  defined  as  the  radius  of  the  vehicle, 
and  it  Is  assumed  that  the  spatial  'correlation  is  essentially  similar  at  launch  for  geo¬ 
metrically  sim'lar  vehicles  and  similar  launch  configurations  when  compared  at  equal 
values  of  k r . 


NoWjthe  naturol  frequency  of  any  mode  of  two  geometrically  slrrrilar  panels  varies 
inversely  with  the  rotio  of  the  panel  dimensions  (or  the  scale  factor).  Consequently, 
since  the  distance  between  nodes  on  the  panel  varies  directly  with  the  scale  factor, 
the  relationship  betv/een  the  scaled  spatial  correlation  of  the  sound  pressure  and  the 
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tcoled  shape  of  the  bending  panel  remains  unaltered  when  compared  on  a  constant 
kr  basis.  It  can  also  be  shown  that  the  geneioli/ed  moss  of  the  scoled  panel  v/ould  be 
the  some  proportion  of  the  total  mass  for  the  same  mode  of  vibration.  Hence,  it  v/ould 
be  expected  that  any  empirical  correlation  for  geometrically  similar  structures  would  be 
a  function  of  kr  ar.J  that  for  this  purpo.se  the  prcvicus  equation  should  be  written 

2  5 

d  u  O  Fff) 

r. _ 

4  W2 


where 


fi  represents  the  constant  of  proportionality 
and  is  assumed  to  be  a  function  of  kr. 


If  is  clear  that  both  the  dampifsg  in  tire  structure  and  the  type  of  structure  will  affect 
any  correlation  of  response  with  forcing  function,  as  both  additional  domping  or 
additional  stiffness  in  a  specific  vehicle  will  reduce  tire  response  for  o  fixed  forcing 
function.  Therefore,  in  absence  of  sufficient  data  describing  these  two  foctors  for  the 
dota  of  various  missiles  in  Reference  3  ,  a  constant  Q  of  15  was  assumed  for  all 

missiles  and  all  natural  frequencies.  The  missile  weight  (W)  wos  taken  os  the  gross 
launch  weight,  but  tlie  stiffness  variable  remains  in  tire  correlotion  as  an  odditional 
unknown  . 

The  forcing  functions  were  calculated  for  surface  launch  of  each  of  the  eight  missiles 
represented  in  Figure  14  using  Figures  8  and  13  from  Reference  3  together  with 
the  appropriate  engine  paromefers,  to  define  the  external  nois"  »*n''ironment .  The 
estimated  value  of  mean  square  force  per  cps  was  obtained  by  determining  the  space 
averag.^.  of  the  predicted  mean  square  pressure  per  cps  on  the  vehicle  and  multiplying 
this  result  by  the  square  of  the  velricle's  surface  area. 
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Thus 


dA 

t 

where  is  the  total  surface  area.  It  is  obvious  that  this  quantity  cannot  represent 
the  true  generalized  force  on  the  vehicle,  l)ut  it  should  be  proportional  to  the  gcner£ 
force  for  constant  values  of  kr. 


A  value  of  ft  was  determined  for  each  of  the  (a/g)  data  points  in  Fif^ure  (14). 

The  results  arc  given  as  a  function  of  kr  in  Figuie  15  .  P  appears  to  be  almost 

constant,  slowly  decicasing  with  increasing  wave  numbet  .  The  standard  deviation 
of  20  log  P  Is  appiox  Imotely  6  db,  v/hich  indicates  that  68  percent  of  the  values 
of  pare  between  0.5  and  2.0,  as  seen  in  the  histogram  of  Figure  16  .  Although  it 
would  be  desiiable  to  effect  a  reduction  in  the  scatter  of  the  values  for  p,  it  is 
actually  somewliot  suipiising  that  tlie  scoftei  is  so  smoll  v/hen  all  the  assumptions 
ore  considered.  Seveiol  comments  me  pertinent,  it'cluding; 

•  First  it  is  improbable  that  all  of  these  missiles  are  dynamically  similar,  since 
a  wide  variety  of  stiuctuial  design  concepts  ore  represented. 

•  Second,  the  forcing  functions  have  been  estimated  through  necessity,  as 
insufficient  measuicd  lounch  noise  data  has  been  made  available.  While  it 
might  be  thought  that  the  use  of  an  pstimator  rather  than  measured  data  would 
help  to  Insure  consistency,  the  estimation  did  not  consider  minor  and  largely 
unknov/n  variations  in  rocket  deflectors  ond  other  configuration  details. 

•  Third,  the  assrrmption  of  constant  C  for  all  missiles  and  all  frequencies  adds 
to  the  deviation  . 


•  Fourth,  it  is  well  !■  rrov/n  that  the  launch  noise  environment  has  a  relatively 
short  duration  and  that  the  response  cannot  be  considered  statiorxrry  from  the 
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statistical  viewpoint.  Consequently,  vaiiotion  in  response  con  be  expected 
at  tlie  same  location  for  diffeient  launches  oi  short  duration  engine  runups. 

•  Fifth,  the  scoffer  at  vorlous  transducer  locations  and  various  flights  of  one 
missile  Is  seen  to  be  almost  as  gieat  as  the  total  scatter  In  tfte  figure  . 


When  these  and  other  factors  are  considered,  the  results  of  this  correlation  ore  en¬ 
couraging  .  Perhaps  the  most  startling  result  Is  the  fact  that  the  meon  value  of 

Is  approximately  one,  as  in  the  single  degree  of  freedom  cose.  This  Indicates 
that  the  rotio  of  force  to  total  vehicle  mass  remains  constant  witli  wave  numijer  and 
that  the  method  of  obtaining  the  forcing  function  used  for  the  correlation  fortui tiously 
gives  the  ratio  exactly  equivalent  to  the  single  degree  of  freedom  case  . 

If  would  be  highly  desirable  to  test  this  empirical  correlation  with  additional  data 
and,  wherever  possible,  to  use  measured  launch  noise  data  In  the  comparison.  If 
would  be  desirable  also  to  test  an  extension  of  these  concepts  to  the  maximum  dynamic 
pressure  portion  of  flight;  however,  the  necessary  flight  parameters  v/erc  not 
available  for  such  an  extension.  In  lierr  of  this  test,  it  oppeais  that  the  value  of  ^ 

In  Figure  16  could  be  used  at  a  desired  confidence  level,  together  with  the  boundary 
layer  pre.ssure  fluctiialion  foi'cin^  function  in  Stretion  .T,  Part  I  of  Refei'cnce  It  to  csLiniato 
this  response. 
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1' N I NK K II INC.  ^ I ( )  1 )  1 ;  K  FOU  KNl  IMHICAF  CORHK  I^\TIOXS 

For  a  linear  struetme,  tlie  mean-scinare  acceleration  res[)onse  at  a  point  can  b(' 
(l('t('rniin('tl  o\;u'f',  ironi  Fcjuations  (l‘jr>)ani]  (120)  for  any  (listril)ntion  of  ranfioni 
proKSuia'S.  In  )r(lei’  to  evaluate'  tliese  expre'saions,  it  is  necessary  to  know  tlu' 
cross-spectral  density  S(  x,  x;  ft  of  tin'  prc'ssurcs  hctu'ci  ii  atiy  two  points  ^ 
wlie're  fluctuating  pre'ssures  are'  applied,  anel  lee  ktiow  the  force-te^-acceleration 
transfe'r  fune'tion  a{\,  y^,;  ft  he'twe'e'n  e*acli  point  x  wheu'e  pressure'  is  a|)|)lie'el  and 
the  pe)int  y^,  whe're'  vileration  re'S[)onse'  is  he'injf  consiek'reel.  Alternatively,  the' 
pressure'  crossee>rre'lal ion  [•(x,  r  )  anel  the'  unit  impulse'  response'  h(x^.  :  r) 

c:ui  1)0  known  and  the  pix'ssiire'  creess-spe'e'tial  ele'nsi'  .nel  the  transfer  function 
can  he  e)l)tain)'d  liy  use'  of  the'  Foiirie'r  transfe)rni  inte'jjjrals  of  Fejuations  (d!))  and  (do). 

Tlu'  transfe'r  fune  tions  are-  le'aeiily  obtained  in  te'rms  e)f  normal  moeles  e)f  the-  structure' 
wht'n  the'  strue'ture'  is  e'le'me'iitarv,  sue'h  as  a  uniform  be'am;  anel  in  such  east's. 

Ke|untion  (IdO)  can  be  re'nelily  e'valurile'd  whe'ii  the  pre'ssure  field  is  kneewt),  eer  can  be' 
dete'r mint'd,  in  terms  of  tlie  pr<'ssuie'  ci-os.s- spectexil  elt'nsilies.  llowe've'r,  foi-  com- 
|)lrx  flight  vehicle  structure's,  the  transb'r  functions  are’  nt't  known  as  ce)ntinuous  func¬ 
tions  of  frceiuency  or  ns  contimu)us  functions  of  position  x  on  the  skin  for  any  other  point 
on  the  structure.  Furthermore,  the  cross -sjx'ctral  densities  of  pressure  are  not 
kneewn  accurately  for  nil  pairs  of  peeints  een  the  skin.  Thus,  in  its  exact  form,  Ee|uation 
(Ida)  is  not  a  practical  engine'ering  te)ol  feer'  preeiicting  the'  vil)ration  respeaisi'  of  a  com- 
ple'x  struc  ture',  nor  is  it  practical  to  u.sc  for  empirical  correlation,  Feju.alion  (125)  can 
be'  use'ful,  howeve'r,  if  apieroxi mations  are  introduced  base'd  on  certain  Icnowledge'  of  the 
prc'ssure  fic'I  I  and  of  the  accci)tance'  and  transmission  characteristics  of  flight  ve'hic'h' 
structures. 
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Vehicle  Surface  Area  Sogmentation 


In  developing  EquoUon  (1?5),  the  spatial  integrations  of  the  pressure  spectral 

densities  and  transfer  functions  were  taken  over  areo  segments  as  shown  in 

Figure  .  A  typical  term  in  this  equation  for  which  m  "  n  involves 

integration  over  o  single  area  segment  A^;  and  this  term  is  equol  to  the 

component  of  mean-squore  occeleiatlon  response  ot  y^  due  to  pressure 

excitation  over  the  area  A  •  If  there  ate  N  such  area  segments/  then  there 

are  N  components  of  mean  squaie  response  at  y^  due  to  excitation  over  each  of 

the  segments;  and  os  seen  in  Equation  the  meon-square  response  components 

2 

ore  summed  linearly.  In  addition  there  ore  N  N  terms  for  which  m  ^  n;  and 

o  typical  teim  is  equal  to  the  crosscorrelation  of  the  two  response  components  at 

^  due  to  excitation  over  the  area  segments  A  and  A  .  The  integrotions 
'  r  ‘  m  n 

were  performed  in  this  manner  in  order  to  account  more  easily  for  pressure  gradients 
over  the  skin,  and  to  take  advantage  of  possible  small  crosscorrelations  of  the 
excitations  over  different  areas  ond  small  crosscorrelations  of  the  response  of 
different  structural  transmission  paths. 

The  significant  soirrees  of  excitation  for  flight  vehicles  arc  the  noise  emanating 
from  the  engine  exhaust  and  aerodynamic  pressure  fluctuations  due  to  boundary 
layer  turbulence,  separated  flow  and  oscillating  shocks.  The  amplitudes  of 
these  pressures  have  significant  gradients  over  the  vehicle  surface;  and  the  regions 
of  high  amplitude  pressures  depend  upon  the  locations  of  the  engines  ond  locations 
of  aerodynamic  discontinuities.  For  example,  in  long  slender  vehicles  having 
engines  at  the  aft  end,  the  sound  pressure  levels  may  be  significantly  lower  at 
the  nose  of  the  vehicle  than  at  the  aft  portions.  For  aircraft  with  wing  mounted 
engines,  the  regions  of  highest  excitations  occur  on  those  portions  of  the  fuselage 
nearest  the  e>  haust  stream.  In  the  case  of  aerospace  vehicles,  localized  excitation 
can  occur  at  the  nose  due  to  turbulence  and  at  aerodynamic  discontinuities,  where 
two  stages  are  joined,  due  to  separated  flow  and  oscillating  shocks.  If  the  sizes 
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of  the  surfoce  area  segments  are  properly  chosen,  the  pressure  level  gradient  can, 
in  some  cases,  l:;e  neglected  over  a  given  area  segment  .  In  other  cases,  such 
os  oscillating  shochs,  the  area  segments  may  be  chosen  as  those  areas  over  which 
the  shocks  produce  significant  pressure  fluctuations. 

ScRiTX -nt  Si7.c  B.is ed  on  Pressu ro  Corre Iji t ion  l,on};th 

Exhaust  noise  end  aerodynamic  pressure  fluctuations  are  known  to  be  random 

piienomena  havitig  limited  sfjotinl  coitelation,  where  tlie  crossconelation  of  the 

pressures  at  two  points  decteoses  ns  the  distance  between  tlie  two  points  increases. 

Thus  if  two  points  oie  sufficiently  seporoted,  the  random  pressures  at  the  two  points 

will  be  almost  sintisticol I y  irrdependent  which  results  irr  effectively  no  crosscorrelation 

and  a  very  small  cirrss -spectral  density.  From  this  starrdpoint,  if  the  areas  are 

chosen  large  enougir,  the-  effective  pressure  correlation  lerrgths  will  be  of  the 

order  of  or  less  tlian  the  dlrr.ensions  c)f  the  area  segments.  With  such  a  choice, 

tire  pressi  res  over  different  area  segments  cs>uld  be  considered  ijncorrelated  so  that 

all  cross«spectral  detrsifles  between  piessurr's  orr  different  area  segments  could  be 

neglected.  Tlie  result  woijld  bf'  tlicM  all  crosscorreinfiotr  teims,  for  which  m  /  n, 

in  Equation  (l/.u'  v/oulrl  l-.e  negligihly  sn'oll;  atrei  the  rrel  response  ot  prrint  y 

would  be  tlie  summation  of  the  mean  sqnstre  response's  due  to  excitation  Over  each 

of  the  area  segments  A  . 

n 

If  tire  surfcjce  area  is  to  be  segrnr'tited  on  tire  basis  of  correlation  It  ngths,  then 
the  si/e  of  the  areas  will  depend  ripon  the  type  of  exritntirrn,  frequency  of  the 
excitation,  and  Irandwidth.  Mer’isured  data  show  that  ocoustic  correlation  lengtfrs 
are  dependent  irpon  wnve-lentjth  itt  the  sense  that  there  is  a  certain  loss  of  cross- 
correlation  of  press'jres  of  two  points  for  eacit  wavelength  of  separation  of  tt._- 
points.  Acoirstic  wave  lengths  at  lav  frequencies  can  be  of  the  order  of  the 
lengtlr  of  a  flight  vehicle,  so  that  in  noria-.v  bandwidths  .rt  lov/  frequencies  the 


aco'jstic  pressure',  are  well  correlated  over  large  surface  areas  of  the  vehicle.  At 
these  frequencies,  the  above  correlation  criteria  for  segmenting  the  surface  might 
lead  to  only  a  few  area  segments.  At  higher  frequencies  however,  tlie  sitoiter 
wove  lengths  are  associated  with  shorter  correlotion  lengths  and  hence,  smaller 
correlation  areas.  Thus  the  higher  tlie  frequency,  the  smaller  these  area  segments 
con  become,  and  the  larger  the  number  of  orea  segments  which  have  zero 

ciosscorielotion  of  pressures.  Correia' 'on  lengths  also  shorten  as  the  bondwidth 
increases  since  the  pressure  signals  tr  nd  more  ‘oword  white  noise  signals.  This 
is  shov/n  in  Appendix  B  v/herc  ihu  autocorr elation  of  fiondwidth  limited  white 
noise  is  discussed.  1  he  deloy  time  r  tiu'  autocor lelotiorr  function  Is  eqr.ral  to 
the  distance  of  sejxrtr ert iroi  ci  tv.'i  piitils  divirled  by  tire  speed  of  propogotlon  of 
the  sound  v.’oves  wherr  the  pnirtts  lii  ul  ing  a  line  jKitalhd  to  the  rllrection  of 
propognti.on  .  In  the  limit,  if  all  ii  •‘quen'^:i  es  vyete  cc>nsideierl  sirnul  taneoirsly, 
it  would  be  found  that  thr’  correlation  lengths  would  bee  'me  vnriobly  small. 
Correlation  lengths  for  '■lerodynnmir  pressure  f lur fuoti '■ns  oic  ver  y  short  compored 
with  acoustic  correlatiorr  lengths,  sin''e  they  ore  of  tire  order  of  the  scole  of  trie 
tirrbrjlent  eddies  In  the  turhulent  lyundory  layer  and  <^et>gtated  flows.  Thrjs  for 
oil  except  the  very  lov.'  frequency  'K<)u'’tir  e.xcitations  the  sri'lnce  rtreo  of  the  vehic 
can  logically  be  subdi  .'i'fed  irrio  ^'orres  o.ru  each  of  which  ill';-  presst,rr'“  rrrnplitndes 
con  be  considerer!  cr'rrstcrnt,  and  bel\v<’en  wlilcb  these  •  hiof  rrriiu)  pr  essrrr ''■s  ron  be 
considered  uncorreigted  so  os  to  hn.e  /ero  volued  cross -spr'ctuil  densities. 

f’c'jrnieril  Sizc  Ba.serl  rm  Structir  r  .'t  I 

Pressure  correlotion  length  anrl  pressure  rirodient  o'e  not  the  only  bases  for 
cboositig  tire  si/es  nf  the  oteg  segments.  Viluuting  energy  frgnsmitteri  olorKj  a  com¬ 
plex  str  uctura  I  [Xi  tb  exp*  r  i*'rtces  I  i  I  t<M  1  try  and  phase  sbl  f  ting  which  '>01  i  es  v/i  tb 
freqrjency;  ctnd  two  structure  transrnlsslon  potbs  within  the  stiuctirre  ran  Irave  qijite 
different  dynamical  cbniacter  istics  .  Random  excitation  supplied  at  tv/o  different 
points  on  the  stirrftirre  vyiil  folio/.'  t /  *>  d i ffet enf  tion'missii m  |intlrs  lit  arriving  at  o 


smcjlo  poitit  wit!iiii  flif  stiijctuie.  If  the  two  poths  Iigvp  quite  different  detailed 
dynurnircil  r hni nctci i stics ,  tiien  the  two  response  components  at  the  single  point 
COM  i  [Kiitinlly  coiielfjted,  oi  nearly  unroi  rein  ted ,  even  thougli  tire  two 
nppli.  1  r  i tctti(Jn^  are  well  correlated.  TItt  dynarriccrl  similarity  of  two  trans- 
Mii  "I  'ti  patle  wirich  slant  to  two  different  points  anci  end  at  a  common  point 
d«  p<  II  Is  in  cienetol  i.ipon  ‘hr  spatial  separation  of  the  two  starting  (xrints,  and  the 
lierp"  -1.  ,•  (inri  hnndwidtlr  ol  the  cycitotion.  Generally,  this  similority  decreases 
•■■I'n  ;n<i*'  rsitig  rlistoncr*  hetweeri  the  starting  points,  '.vith  Incteasincj  frequency 
and  v;il!i  mrresising  hrjndwid*''.  Irr  this  Sense,  the  ctossroi  relations  of  responses 
rtlonn  two  struf  final  tr '•insmi  ssisrn  paths  have  charar  ter  isfi  rs  v/hich  ore  similar  to 
thr*  f  r  nsss  frr  r  #*  1- 1  r  t '  I  MS  f  ,it  trie  'ipydtea  tan'iiim  ptr'ssurr's. 

Ih'pitliir  Hi  l■!l(■(lIl'''s  l  enpatul  pfitnafilv  in  wr'll  dt-fir.i'd  ini)t|<>s  nf  vilu-nl  ir  m,  .-mtl 
the  enri'el  il  iiMis  of  icsp, nrse.s  ;t(  .r  point  diir*  to  irrpnl;'  .tl  l\so  ililTi  i-ctil  points  tn:i'. 
he  well  eor  t-ehtlr  'I  if  the  inprrts  :rrr>  corTr  intr  d.  In  n  eonip!e\  si iniclni’e,  the 
rr'sponse  (ml  hs  ,t  ne  nittch  inoi-e  eonqrlex;  Itowr'vr'f,  the  lack  of  erosseoirokrlion 
hetwer'tr  resirintsi  s  alonp  two  paths  repn'sr'itl s  a  simiiticant  sinrirlifiealion. 

'I'htis  the  siiffaee  at'ea  sejrmr'nts  can  ai.sn  Ire  chosen  en  the  basis  of  si nictiinal 
r'oi'i'r'lation  h’lrnlhs  srt  tfrat  Ihr-  rTosscorrclat ion  r'‘.s|>onse  teians  in  lirinalion 
(ljr»)  can  hr-  ar'ph'clcil.  tt  is  to  !»<•  noted  that  the  si  i  iiclii tar  I  crri  ri'lat  ion  Icnptlis 
may  Ire  slntrtrM'  thc.ir  tlrr-  (rressurc  cora’clation  lengths,  such  .as  in  th<‘  case  rrf 
a  stiflencil  cylitrdef  rrr  .stiffenerf  (ranel  (wciti’d  hv  acrriislic  wavr'S  whcaa'  th<' 
(lislatrcr's  Irelwcrnr  sliffr  tiers  is  much  less  lhait  lire  \'av<’  Ir-npllr  rtf  the  sttutrti 
waves. 
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Response  Spectral  Deniuty  with  Crosscorrolation  Terms  Ncj^lectt  <1 


Based  on  ttte  above  comments,  it  is  reasonable  to  assume  that  Equation  (1?r>)  for 
tt’o  response  spectral  density  cari  be  written  in  the  simpler  forrts 


S  (f) 


r>i 


Z//  aV,  a  S(?,?  0  dA(?)  dA(j’) 


n  1  A  A 
n  n 


(1:55) 


The  two  integrations  In  fqnatior''  (  )  '5'  are  taken  ovru  a  common  area  segment 
A  shov.n  in  Figure  (11)  ;  and  tire  resulting  term  in  Equation  (1^.51  is  equal  to 

tire  aereletatlf)n  resporViC  spectiat  ^!•'rt^l  t  y  ot  due  to  flucttjoiinq  pressure 

r 

('xcitatlnn  rjver  the  area  A  .  lire  net  spectrol  density  of  rc'ponse  at  is 

n  ■ 

then  tire  sumrt^atIf)n  srjc.lr  response  spectral  density  components  over  all  N  area 

segnrertts  . 


Before  irroceedlng ,  '  •  should  be  noted  tl.nt  for  o  given  frequency  and  [bandwidth, 
the  area  segments  chosen  will  In  general  be  diffetfuit  for  different  types  of 
excitatiorr;  and  hence,  'sgxrrctfe  on^rlysi’s  shotjid  be  conducted  for  exharrst  noise 
and  ftsr  or^rodynarriic  frressure  f luctrrations  . 


It  is  crssumerl  now  tii^r’  the  atr'n  seqrrrf'rrts  have  been  so  r.  ‘Osen  that  the  ptesstjre 
spec  tral  density  at  n  qi  .'fn  freqrj<w,r/  i-,  ronstunt  over  *'och  area  seqrnent  .  De¬ 
noting  the  pressure  spectral  dr-nslty  over  A  as  S  (ft,  Eqtjollon  f  Idoi  con  be 

n  n 

written  irr  the  form 


S  ff) 


N 

i:  h."’ 

n  1 


ff 

A  A 
n  n 


♦  -» 


f)  a  (> 


— f  d  — r ' 

y  ;t)  C(x,x  ; 


fV)  Afy  )dA(x 


wliere  Cfx,  x  ;f'  iso  spatial  correl.rliorr  coef f i ci t'nl  of  tire  applied  pressures 
which  has  rr  rnoxirrr  m,  (jlsolutf>  voliro  of  irnity.  T  hi  s  coi  t  e  lot  ion  coef  f  i  ci  eri  t  whicl'i 


^7 


is  defined  b/  the  equation 


mJIk  * 

S(7,K;f) 

sin 

n 


(137) 


is  equal  to  unity  whore  x  x  and  is  approximately  equal  to  zero  when 
I  -  "x*'  j  is  large.  A  typical  voriation  of  c!?,?  ;f)  v/ith  separation  distance 
x‘|  fot  a  homogeneous  pr<‘ssuie  field  is  given  by  the  exponentially 


X 


damped  cosine  function, 


•d(f) 


cos  k  (f)  •  I  X  -  X 


(13S) 


where  tlie  constonts  d(f)  and  k(f')  aie  fif'quency  dependent .  For  a  plone 
wove  acoustic  field  at  fieqnency  f,  t!\e  constant  d(f)  Is  equol  to  zero  and  k  (f )  ■ 
wheie  c  is  the  velot.  ity  of  propagotion  of  tlie  waves  in  tlic  direction  of  the  vect 

X  -  "x^  . 

The  double  area  integral  In  Equotlon  (  l.■5."))  is  of  the  form  of  the  square  of  an 
average  force-to-occeleiation  transfer  function  between  the  area  A  and  the 

point  y  .  This  can  be  seen  by  notitig  that  S  (ft  dAfxt  dA(x')  has  the  urrits  of 

2\  2 

(force)  /  cps,  S^(f)  has  the  units  of  (acceleration)  /cps,  and  by  comparing 

Equation  (l.T."))  with  the  definition  of  the  transfer  function  given  it)  Equation  ( 

Thus  if  a  denotes  this  aveiarje  transfer  function  at  frequency  f,  then  tf 

acceleration  response  spectral  density  hr  somes 


N 


S  (f) 


L 


A  ^  a  ^(^;f)  S  (f)  (133) 

n  n  r  n 
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'.•/fiorc 


V 


a  7y;f) 
n  r 


f  j  ^  (^'/y'^;0  C(x,x';f)  dA(x)  dA(x')  (I  ***) 

•'a  -/a 


n  7^  TV 
n  n 


'I'lie  Avcrat^c  Transl'ci  rimotion  in  Tenns  jf  Frc-nioncy  B:tn(}s 

f  ni  ('Irniontai  y  stnicfiH  ('■>  linviny  pinsMiio  whicli  coti  he  defined  in  teims  of  the 

coiielation  coefficient  C(x,  x';f),  the  oveioge  tionsfei  function 
'•tr'tf'i nii ned  in  teims  of  nofinol  tnode«.  of  the  stiiictuie.  For  coupled  structure 
this  '■j'/einge  tmnsfei  func  tion  cntronly  he  determined  r'xpetinrcntally  .  Average 
liomtf'i  f'.m.ctions  of  this  lyfre  vjeio  rteteimined  expei  imeritol  ly  for  the  Snot!' 
missile  os  ie|)Oilerj  in  Reference  (IPt  .  Measured  \ihint!nn  dotrj  are  trsirnlly 
token  in  freeprerrey  honds,  sudt  os  1  3  -  octave  limids,  ctnd  it  is  of  interest 
to  'fevelop  on  expiession  for  the  meon-squcrrc;  ncceleiollon  tespemse  in  such 


lir’quency  Ironrls  ,  let  V<  fy;f  ,Af)  drurote  tie. 


•  meon 


-squore  acceleration  re- 


'ponse  at  the  point  y  in  tlu'  fiequetrcy  bond  At  cr'ntei.-d  on  lire  frequency  f  . 

I  c 

lire  meaci -sqiKir  e  respi'inse  Is  giverr  by  tire  equotlf'ir 

f 


••2  ‘ 

v;  (y  ;f  ,Afl 

r  c 


f 


S  ff) 

t 


df 


(1 


wlif'te  f^  and  f  denotes  the  lower  ond  rrpper  freqiipnc,  limit'-  ol  the  hnnri 

(Af  f  -f ,1.  If  the  pressirre  'pectrol  density  S  (f)  is  tea'-onahly  corrstant 
ir  L  '  n 

over  the  frequency  iKitrd  Af,  and  if  S  fl  ,Af)  clertot"  .  rl.i'  overcroe  value  of 


n  c 


S  ^f  )  in  ihis  frarirl,  then  suhstitut loll  of  Krpialioii  (  in!))  into  (Ml)  ^;ives 
n 


V/  (y  ;f  ,/,f! 
r  c 


N  - ; - 

fAf)  V  A  ■  o  (y  ;f  ,AA  .  I  (f  ,Af)  (I  f:!\ 
/  j  n  n  '  I  c  n  c 

rt  1 


wlteie  »Ke  quantity  a  Af)  is  the  space  overage  and  frequency  average 

value  of  the  tinnsfei  function  he^//ee^  the  orea  A  and  the  point  "y  ;  and  this 

n  r 

quantity  is  defined  by  the  equation 


a  Iy;f  ,A0 
n  'r  c 


(n:i) 


K(lii:itian  (112)  h:is  :»  niorr'  useful  for  in  for  practical  engineering  application  than 

its  exact  cnunleipail  given  by  the  (> spression  iinmecliately  following  Efiuation 

(120).  The  only  assumptions  that  have  been  made  in  tbveloping  Eriuation  (142) 

are  that  tlio  pressurt'  spectral  deasities  are  constant  rrver  each  of  the  surface 

area  segments  A  ,  and  that  the  s'/t's  of  the  area  segments  were  so  chosen  (foi 
n 

the  fr('f|uencv  band  Af  centeia'd  on  I  rmiuency  f^_  )  that  (’ither  pressures  acting 
on  the  surface'  are  ancorrclated  ln'tween  <liffcrent  segments  and  A^^  (m  /  n) 

or  the  responses  at  T^.  d  ie  to  exciiations  over  diffro-enl  area  segnu'nis  are 
uneorrehPed. 

L'  y  h  !  it  i-'i* A.Vil*li.‘ii  •*.  Xl'i.L'b  I  ‘^’.11 T’, '  L15-4 ! 

At  the  piesor\t  state -.<1 -the -oi  t .  tiov,evci  .  it  i*-  nerrs'.oiy  to  oppiox  i  mat*-  the 

aveiriae  tnnsfei  fotictian  a  (y  ;f  ,Afi.  To  shoiten  the  notation  let  this 

n  I  c 

average  ftntisfci  fuc.ctlon  lx  denoted  !)y  the  symlio I  a  .  In  oide*  to  niiive  nt 

a  leasofinhle  englof'eiiog  model  foi  reptesenting  a^  the  stiuctuial  tiansmission 

path  vdll  he  assumed  to  he  mode  up  of  three  series  of  elemerrts  as  sliown  in 

ficjuie  (12)  .  Tlie  (i;  t  rdeni'a  l  i^'piesents  the-  acceptonre  hy  the  stiuctuie  rrf  the 

?  . 

incidence  suifcrce  pif'S'-ni'-s  om)  it  rj  mean  squnie  force  1  (l),  hoving  n  constant 

n 

spectral  density,  Is  applied  over  the  area  A  ,  tlien  the  mean-square  acceleration 
response  of  fire  surfcrce  str  rclute  vyill  he  denoted  by  tbe  quantity 

a  ,Af)  •  f  ^ft) 

n  r  n 
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wlicie  3  is  o  spcce  ovcroge  and  froqucncy  overage  tronsfer  function  of  the 

2  . 

suifoce  structure  over  A  •  The  quontity  B  contains  the  effect  of  joint- 

n  n 

acceptance  of  skin  panels  at  the  surface  of  the  vehicle.  The  second  element 
in  the  series  denotes  the  transmission  of  this  vibrator'/  motion  from  the  area  f' 

r 

to  the  region  of  the  point  ;  ond  the  output  of  this  element  is  denoted  by 


T 


>  ^ff  ,A()-  3  ,An.  F  ‘(If 
nr  c  n  c  n 


TItrjs  the  a/ celerntlon -to -'^tcceh'rot’on  transfer  function  ")  gives  the  rotio  of 

nr 

tire  mean-square  response  'Ift  the  ttoiqhlsorhood "  of  y  to  the  meon-squate 
acceleration  response  ''in-tli<?-neigiihoihood'*of  A  •  Experiments  on  the 
fuselage  of  the  Snath  missile,  ir 


(vsi  ted  in  Reference  (12)  shov.  that  q 


nr 


con  be  ftp[)t  ox  I  mated  by  /in  »■' 


itiol  decay  function  with  dis*ancc  olong  the 


structure  as  the  atgumerrt  of  the  function.  Thus 


V  ^(f  ,Afl  e 

nr  c 


8(f  ,Afi 


(Iff) 


un 


wlrere  1/  in  this  equation  denotes  a  rxsition  of  a  tepres'‘nt<'tl  r'  jv'int  withi 
n 

the  area  A  .  The  qncjntity  8(f  ,Af'  is  a  constant  whi.  ii  dj^penrk  up/rri  bnndwidtli 
n  c 

atrf!  frequency.  Tire  fin/rl  elemen*  showrr  in  Figure  (!.')  drurotes  tlic  or.'  eptance  of 
strijctmally  tr ansmi tte d  ener /g /  by  a  "receiving''  element  whi  h  may  be  n  structural 
component  or  a  piece  of  equipment,  letting  tfr  (f,  Af  )  dcfiote  tire  accelerotion- 
to-accelercitlon  transfer  fiin/  tion  at  the  "receiver  the  re-’on-squaie  acceleration 
response  of  tlie  receiver  is 


tJi  ^(f  ;  Aft  •  d  ^li  ,  At)  •  3  ^(f  ,  Af)  •  F  ^  (t) 
,  c  iM  ^  c  n 


It  follows  that  the  average  transfer  function  a  ^(y^!  ^0  can  be  approximated 

by  the  product 


a  (y  ;  f  ,  Af)  5  P  '(f  ,  Af) 
n  r  c  n  c 


■»  ^(f  ,  Af)  •  »  ^(f  ,Af)  <’  *■’» 
nr  c  r  c 


(f^,  Af) 


space  average,  frequency  average  force-to-acceleratior 
acceptance  function 


y 

nr 


^fc.  Af) 


fiequency  aveiuye,  occeleiullon-lo-acceleration 
transmission  function 


4.  ,  Af) 

n  c 


fiequency  average  acceleiation-to-occeleration 
receive!  function 


where  AO  is  siven  niiproxlinately  by  Equation  (Ml), 


Suiistitutinj.',  Equations  (I  15)  ii.to  (112)  j;iv('S  the  final  approximate  expression  for 
the  mcan-sriua'-e  response  in  a  <;iven  freqneney  hand  in  terms  of  the  aei-eplance. 
t ransmission  and  receiver  functions: 


n  I 


where 


w,  -  VV^  (y_.;  t^,  A|.) 


?  (f  , 
n'  c 


Aj.) 


2 

The  quantity  is  the  approximate  mean-sriuare  force  acting  on  the  area  segment 
.  All  fjuantities  in  Equation  (1  IG)  arc  averaged  over  t1\e  frequency  Ixind  /J  whii 
has  a  center  fre-iuency  f^  .  'I'hc  iivlox  n  refers  to  apiilied  excitation  on  the  n-th 
cirmpartnienl. 
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Dctcrminat i on  of  Tran sni i ssion,  Receiver,  and  Acceptance 


The  seoarotion  of  the  transmission  function  >  and  the  "receiver"  function 
2  nr 

i|)^  is  in  order  to  account  for  local  structural  parameters,  sucji  as  local  inass 

and  stiffness,  of  the  "receiver".  In  this  connection  a  us'^ful  concept  is  obtained 

b/  dividing  the  structure  into  bays,  as  shown  In  Figure  (18),  letting  > 

account  for  tlie  average  axial  decay  in  response  level  between  the  region  A 

wliere  excitation  is  applied  and  the  boy  which  contains  the  point  "y  ,  and 

2  2  2  ^ 

denoting  •  F  (t)  as  an  average  response  of  tl.e  structure  in  the  bay 

in  wlilch  ~y  is  located.  The  product  of  this  average  bay  response  and  the  receiver 

.  'T  . 

function  tji^  then  gives  the  response  of  the  particular  structural  or  equipment 
component  located  at  .  Also  shown  in  Figure  (18)is  a  set  of  acceleration 
response  points  for  various  structural  components  along  the  axis  of  the  vehicle, 
and  it  is  seen  that  the  response  dota  are  scattered  about  tlie  exponential  decay. 

The  differerence  between  the  exponential  decay  curve  and  each  of  these  response 
points  is  accounted  for  fry  the  receiver  function 


2  2  2 

Experimentally ,  the  functions  B  ,  7  ,  4>  can  be  determined  by  providing 

n  nr  n 

localized  acoustic  excitation  to  tli"  vehicle  skin  in  a  metnner  similar  to  that  done 


for  the  Sndrk  in  Reference  (12)  .  for  convenience,  the  vehicle  can  be  segmented 

along  the  axis  Irito  strurturr:!  zoun.;  otrei  for  compai tmented  vehicles  leaving 
bulkheads  or  other  material  structural  divisions,  these  zoires  can  bp  clmscn  a.s 
the  compartments.  Acoustic  excitation  could  be  applied  to  a  portion  of  the  skin 
of  o  single  compartment  ctnd  ncceleratlor.  response  meosurements  made  at  a  number 
of  points  witidn  the  same  rompni  f  ri' '■nt  and  in  each  of  the  other  comportments.  By 
averaging  the  mean -squtrt e  acceleration  responses  within  each  of  these  compart¬ 
ments,  an  average  resoonse  of  eacli  core (Xirtment  could  be  obtcrlned.  The  larger 
the  number  of  points  within  the  compr!r*ment  ond  the  wider  the  range  of  structural 
compcjnents  consideed,  the  mute  representative  this  overage  response  will  be. 


When  tlicse  average  compartment  responses,  in  a  given  frequency  band,  ate  plotted 
against  oxini  distance  along  the  veliicle,  it  will  Ire  found  that  the  aveioge  accel¬ 
eration  ampliliicles  will  decrease  with  distcrnce  from  the  comportment  directly 
excited  by  sound.  By  fitting  an  exponential  curve  of  tlie  form  given  by  Equation 
(bU)  ,  to  tliese  data  points,  the  parameter  6  (f  Af)  in  Equation  (bt4)  can 

be  determined,  and  from  tiris  the  trarv-missron  function  7  will  be  knov/n 

nr 

from  liie  directly  excited  bay  to  any  other  bay.  Reporting  this  process  by  directly 

excitlrtg  (.-uch  irf  tire  Individual  compartments  and  determining  average  responses 

- - 

for  '■ill  comportmmrt'. ,  the  complete  sr-t  of  >  (f  i  Afffunctians  between  the  n-tli  and 

nr  c 

and  r -til  compartments  can  be  determineri.  Clearly  by  use  of  filtered  response 

2 

data,  the  firnctlat>«,  )  ran  be  obtained  for  each  frequency  band  of  Interest. 


It  has  lieen  lucitly  rrssnmed  tliat  tire  hinclions  t  are  equal  to  unity;  and  in 
tills  sense,  •>  represents  the  axial  rler<iy  of  overage  response  amplittx4e  between 


two  different  comfxrr  tments  . 


An  average  acceptance  function  p_  (f  ,  Af)  can  also  pe  obtained,  for  the 

n  c  2 

n-th  compartment,  from  those  meosuted  data,  since  S  can  be  defined  os  the 

n 

rot:''  of  the  avf'toe.e  mean-sq'Kite  oc- elerotion  response  in  tlie  compartment  to 

tlie  totol  mean-square  exci  tatlon  \  octing  on  llu'  suil.ret.'  iitc'.t  uf  thul 

2 

compartment.  He  e  also  P  cart  Ire  determined  for  eoch  freqrrency  ban  !  of 

interest,  for  most  velticles.  It  con  he  ossimierJ  tliot  If  the  excitation  is  applied  to 

only  a  portion  of  the  surface  area  of  n  given  compartment,  then^if  this  portion 

of  sutffree  structure  directly  excited  is  t ent osento ti ve  of  the  entire  sijrfnre  structuve 

2 

of  the  compnt trnent ,  the  acceptance  funrlion  P  will  he  representative  of  that 
obtained  for  ony  other  portion  or  the  whole  of  the  surfoce  of  the  com[iortmcnt . 


2 

It  should  he  tutrerl  that  th*'  ohovr*  efefinitions  of  the  acceptance  function  p 

and  the  transmission  function  >  differ  shglitly  from  tlie  oricjlnol  definitions 

nr  ' 

v.liicli  hni-ert  these  functions  directly  *•'  surtoce  strrrcture  ratlicr  them  to  the 
str  ijr  tiit'-il  7orie  <>i  r,c)mpr-tr  tro'iit .  At  the  present  stote-of-the-rirt,  however,  these 
nt»  !  ortsorKihle  llrst  npptox  linrtions  r-,f  tlirii  tnoie  exact  counterparts.  As  rid  litloncrl 
tmpr  'Iorif  >'  drrto  on-l  other  reensirred  vibration  response  doto  ore  obtoined,  there 
I'uneliutis  e.in  nritu t  rtl l\  he  uivr’n  more  ri  fitrc  rl  trentnumts. 


94 


The  Si^nificiinco  of  Loc  il  ^lti'vicUn':il  Paiumclcry! 

Ilic  rliflf  i»,‘iKe  liftwecn  llie  mecjn -square  respcinsc  of  o  structurcil  element  cind  the 
uvriuq'’  mron-wjDOie  le'.ivjnse  ot  the  ho.'  In  v«hlch  the  eicnictit  is  locfited  is  equal 
to  the  receiver  function  U'  d  '  for  that  element.  This  function,  which 

I  r. 

depenri'.  (iiimfiiilv  upon  ioi  (j|  tiioss,  stiffness  and  dompiny  prerperties  of  the  elcrmcnt, 
shorrid  he  app'ovim'iiel  /  the  srieie  lo:  cfiiert  excl taticrii  of  tlie  comportment  con- 
fui'inii  iIk'  element  nnU  e/t  i  trrtlori  of  any  other  comprjt  tment .  At  present,  there 
i s  'ei  >  !  i  tt  i  ‘  i-  tr  ■) wn  1,0:1''  ei  ni nq  l!  dislr  if  ution  of  etioiq/  thr  ouqhout  a  complex 
strircluii  ,  pur  ticulor  ly  wi'l.in  o  lorrilired  tegiorr  siicli  os  o  comportment  where 
I'.rul  siiu'*Mr'il  d(  lo  I  si"  iifica'itl ,  qileet  res|)on5r<'  ninpHtudes.  In  ofdef  lo  sftidy 
this  rispr'f:t  '  f  tlo'  strur  turnl  y'ht'jti'''n  pi'i!)lem.  Statistical  onolyses  can  lie  mocle 
in  whi.-h  the  tt'ceivet  function  (Ji  is  correlated  v/ith  such  local  stttjctjrol  parametr'r 

as  moss,  sfi ffness,  lengtii ,  choroctei  istic  or  effective  thickness  for  flexurol  vibrations 
of  structvjtal  component,  rotio  of  moss  of  component  to  moss  of  neighboring  components, 
(similorly  for  stiffness  rritlo'-t,  curvature  '>f  element  as  in  the  case  of  skin  (Xinels, 
or  i''ntcttion  relative  tu  tf  <■  primary  <fIi<'f.tion  r>l  comportrient  r'ycitntion,  ftequetiC), 
liondwifltir ,  atttf  trumlrer  of  resporises  In  a  fiequr-nry  bard.  Other  significant  local 
parameters  tnoy  also  become  evident  r'utltr'j  such  a  stotlstirol  strrdy  .  It  is  felt  that 
a  knowledge  of  such  cor  relatioris  worrld  be  of  significant  value  in  the  development 
of  an  eerplriccilly  hosed  vllurttion  te'prrttse  ptr'dictifui  technlqrre  . 

By  obtaining  nrea'vrred  vilurtlion  datcr  on  a  numl'et  of  fllgfit  vel'Icles  of  tlie  sorr’e 
'-las,  and  of  n  number  of  cflffeient  classes  of  vehicles,  statistical  onol yst^s  could 
he  corrdurted  to  stu'^h,  the  accet'taiii  e  ttnd  ft ntts'rt'ssion  functions  R 

n  nr 

Thr-  tirnrtlons  cotrl'l  fse  C'3i t ela ted  witfr  m'^te  gt'oss  vefiicle  parameters  such  as 
nrass ,  lencjfli,  torfitrs  crnd  itveroje  surfcjce  density  of  a  compartment,  'ivetall 
i»tii.ttlt  rrrrd  vyei'iht  uf  a  vehicle,  f  teqr  re  iic /,  I  anrlwifltli ,  a/eta<je  drimpln'.'j, 
stxirirrl  ''iii»>|  iticui  leuaths  f'-r  S'jrla  r-  [>resstjtes,  per  emit  st' artirrcrl  weierht  of 
t'ltol  vehicle  welgfrt,  etc. 


Much  could  be  learned  from  iuch  statistical  analyses  of  complex  structuies 
including  the  'caling  lows  which  govern  the  response  levels  of  diffeient  si/o  but 
geometrically  and  dynamically  similar  vehicle  structures.  Scaling  laws  for 
exactly  similar  structures  con  be  obtained  easily  from  fundamental  dynamics 
concepts;  however ,  hecnuse  of  the  interplay  of  the  host  of  overall  and  local 
vehicle  parameteis,  the  exact  elementary  scaling  laws  provide  only  a  first  order 
approximcjlion  of  the  true  scaling  laws  for  complex  structures. 

l.inHt:iUon.s  on  tiro  Tsc  of  DrtJ^.i 

In  Ronoi'.'il,  nio.tsm-cd  cNcilrttion  .ond  iv.^ponsc  rlol.t  fronr  flij’hls  oaruiol  Ire  used  diree 
to  oldoin  the  acceptance,  transmission  ond  receiver  functions  in  Equotiorr  (1  Id) 
since  there  is  insufficient  I'Inta  for  flie  rtnmhet  of  unfnov/n  constonts  to  be  determined. 

To  show  this,  let  the  suhsciipt  m  denote  points  v/itliln  the  fn-th  comportment  at 
wirich  tesponsr>s  hove  iieen  mecisuied.  TIren  tire  response  ot  tire  r -th  point  in  tire 
m-th  compartment  Is  cjlven  b/  tin*  *'xpiession 


Tire  Ofceptonce  funct!or>  is  assr.r  in  ferl  with  excitation  o.in  lh('  n  'ii  omprntment; 

is  the  fractlonnl  loss  oi  irieon-snuaie  ar(ehMnlIon  lesp-  n.  ■  between  the  n-th 
nm  2 

and  m-th  comportments;  ilj  is  the  leceivet  function  for  the  r  -  th  [-.olnt  In  the  m  -th 
— y-  r  m  ' 

comportment;  a  .d  is  the  mean-square  foire  octinc)  on  the  m-th  r  ■jmp..it  tment  . 


where 


A  X  -  X  overage  distance  betv/eeri  compartments 

n  m  n  m 

8  8(f  ,A^)  "  constant  to  be  determined  . 


TbeiP  is  one  such  equation  for  each  compat  tment  (m  "  1,2,3/...,  N)  so  tliot  thei  e 
aie  N  sucli  equation'  .  There  are  N  values  of  N  values  of  and  one 

value  of  8;  and  hence  tliete  aie  2N  '  1  unknowns  to  be  detormitred  f.om  N  equations. 
Tims,  even  when  the  'icceptonre  and  tirmsmission  functiorrs  are  referred  to  compartments, 
anri  the  ttatrsmi ssion  function  is  appio-. imoted  l>y  tire  one  parameter  exponential 
function,  there  is  insufflr  lent  dnto  to  determine  the  unknown  quantities. 

It  is  possible  to  obtain  'it  least  *iisl  older  appio'^imotloris  of  the  unknowns  appearing 
in  Equation  (I.ul),  |t  'r.ensuie'i  cl  lotiopi  and  vibiotlon  response  data  are  ovollable 
for  N  *  1  clltteient  flirpa  (  or  qi'mcid  run  --  up  or  groui'it  firincjl  r.orrdi  t  ioirs  .  V/hon  the 

distribution  of  oonlled  pressures  over  tire  vehicle  Si.)ifri:e  chanqes,  then  the  menn- 

■)  .  > 
siiirnte  f'.icr  s  f'“  \:ii  \,  :r  i|  st,  d  tli<  .vr >  rcc  :i v 'r'.'ipu  I'esp'  isi  :-  \''i  .  Iiii.'-',  it  i 

I  he  N  1  (|  ( flV  re.it  fl  t  e  .1(1  ii  i  ,rs  I  l,ei'('  ;i  ic  N  J  d  it  |e  i  i-  it  lu  '  'SS'i  i  >  ’  d  i  si  i  i  I  nit  i  , 

Over  the  sur^oce,  tlien  there  nts’  N  •  1  equations  of  the  for  m  ( 1  ■>  i )  for  euu.li  vaim;  m, 

~T  .. 

and  hence  there  ore  ?N  •  1  equations  in  th.e  unk nov/ps  'T'  ,  P  ,  h  I*  is  o'sulne(^ 

-y  rtt  ri 

of  course  that  tlu'  avera')''d  accerrtnnce  furirtlons  P  are  unchantred  fer  dir.  diftr'ent 

n 

force  distr  lljutloris  .  This  appears  to  Ise  a  lensnnnble  n[ipi  ox  I  nrrj  I !  rr  for  o  cp/en  type 
erf  evritotioo  In  which  he  pre'S’.irf'  sprite  cor  relri  tlons  at  the  '•ehj_;._l/’  surfure  lerroin 
oppro-ximatel  V  fixed  ,  Otice  numerical  values  of  the  qi'antitles  p  and  8  cic'  i  no'.'o, 

i  .  ' ' 

the  overoqed  leceisfu  functions  iir  can  Ire  determined  h.)i  r  cich  s ||  ijr  t  y  ol  ■  rnc[->one’', i 

1  m 

hy  use  cjf  Equation  (1  \  . 


A  final  slmpllf  Icai'on  uf  Equaticni  (1  l<'>)  i,  now  tonsisln  s’d  .  The  purduct  (■'■  •  [ 

an  approximation  of  the  rj  -f'lriqe  rerron -sq- roi e  acceierotion  lespio  ise  c,  f  the  n-th 


compartment,  and  as  such,  this  product  should  Kt-  inversely  proportional  to  the  square 

2  ~~7~ 

of  tlic  total  mass  M  of  that  compartment.  Thus  the  quantity  (M  P  )  should  be 
n  '  ^  •  n  n 

nearly  independent  of  the  total  mass  of  the  compartment.  Fuitheimore,  if  the  surface 
structures  of  all,  or  most,  of  the  vehicle  compartments  are  similar,  the  average 
joint-acceptance  (or  acoustic  stiuctural  coupling  between  the  applied  pressured 
and  the  veliicle  shin)  sliould  be  approximately  the  same  foi  eoch  compartment. 


Since  the  acceptance  function  p  is  clearly  proportional  to  joint-acceptance,  tlie 

2  T"  .  ” 

quantity  M  p  micpit  have  opproximatel y  t!ie  same  value  for  all  compoi  rments  . 
V^ith  this  assumption,  Fquotion  (l.ul)  aon  lie  written  In  the  form 


b  —  . 

n  n 


(ir)2) 


w'here  D  is  a  constont  to  he  cieternuried  olonq  witli  tlie  quon’itirrs  rj-  ond  F.  Therf 

nr 

ate  nov/  N  *  2  unhnovvns,  so  tliot  in  addition  to  n  ensured  data  irr  eoch  of  the  N 

cr’^mpartments ,  it  Is  necr'ssoty  to  have  n!<*a'tiied  data  in  two  r '.'mpor ‘'nruits  for  o  seconf 

flight  condition  hovlnfi  a  diffei‘'“nt  iti'.tr  ilnrtion  of  fluctuatlno  pressures.  ()nne  D  ont 

2 

8  are  known,  the  receiver  functioru.  i|)  con  he  founrl  from  the  eq  lotion 


-A  f> 

o  r  '  M 

II  n 


|ri  conclusiori,  Equotions  a-ul  (ir>l)  (lali)  mn  be  ised  \n  cieteimiiie  e'noi r  1  r'l 1 1  . 

occeptonce,  transmission  ar.rj  local  vibratory  cfraroc  ter  isti  cs  of  rompley  struct, j:e' 
if  thern  is  sufficient  niepMitorf  ciftta  t(j  ffesetibe  tire  r'^ritation  ond  the  resportse.  Tire 
veliicle  con  be  sc^gmentecj  Inro  cornpot  tments  os  dc-'-tihed  above,  and  it  iv  necessary 


to  df  'erminp  the  average  meorr  square  force  acting  on  the  exterior  of  each  cornpoit- 
menf  and  to  have  ovollabte  a  representative  sarnpling  of  vibration  response  data 
witliln  eaclr  compartment.  It  is  nccessnr/  tlrot  sucli  measured  data  be  available 
for  a  number  of  different  flight  conditions  in  which  the  mean-squoto  forces  over 
the  segmented  areas  of  tire  vehicle  have  different  distributions.  V.'hcn  data  of  this 
type  is  not  nvallahlc,  flight  data  cannot  used  to  develop  empirical  corrclatloi 
ptedlf  tioir  methods;  ernct  in  tlris  case  it  is  necessary  that  imperlonce  measurements 
tor  firry  equivalent  mecisurertrentsl  be  user!  to  evaluate  tire  acceptance,  transmission, 
and  lof;al  d/fiamic  characteristics  of  vehicle  structures. 
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APPENDIX  A 


AUTOCORRELATION  FUNaiON  FOR  V/HITE  NOISE 


Let  x(t)  be  a  signal  resulting  from  a  stationary  and  ergodic  rondom  process  which 
generates  white  noise.  The  autocotrelation  function,  I  ),  of  x(t)  is  defined  as 


r(i  )  -  Lirn 

T  •  00 


]_ 

T 


(t)  •  x(t  ft  )  •  dt 


(A.l) 


Vv'hen  x(t)  represents  white  noise  containing  all  frequencies,  the  autocorrelation 
function  rfi  )  is  proportlotxil  to  the  Dliac  Delta  function,  5(i  ),  which  is  defined  as 


S(t  )  00,  t  0 

0 , 1  /  0 


f 

0 


6(t  )  •  di 


-  f 


6(t  )  •  di 


1,  2  ♦  0 


(A.  2^ 


foi  any  f  >0  no  mattei  how  small.  Tliis  fact  is  not  dliecfly  obvious  fiom  (A.  1).  I  i 
oidei  to  prove  this  point,  and  in  order  to  determine  the  constant  of  piopoi  tlofial  ity , 
it  is  necessary  to  us'^  the  constant  power  spectral  density  definition  of  wliite  noise 
ond  the  V/lenet  -  Khlnchln  relations. 


Let  Sfw)  denote  the  power  spectral  density  of  the  signal  x(t).  Since  x(t)  is  characterized 
by  the  properties  of  v.'hite  noise,  the  function  S{u)  constant.  From  page  67  of 
Refeietice  (1A)  ,  the  '.Vlenet -Kiilnchin  relations  are 
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(A. 3) 


(A.  4) 


The  symmetric,  or  cosine,  forms  of  these  relations  ore  used  because  the  functions  S((j) 
end  I  fi  )  ore  symmetric,  tliot  is 


l(-r)  I(t) 

SM  S(u) 

The  latter  relations  are  given  Ly  Equation  (1-61),poge  20,  and  Equation  (2-1  5), page  3 
of  Reference  (1 01 . 


At  the  upper  limit,  the  function  In  (A.5i'is  indeterminate  since  the  sinu.T  can  lie 
anywhere  in  the  range  -1  <  sinW"^^  1.  Neglecting  this  indeterminate  upfX"i-  limit, 
r  (t  )  is  seen  to  eriual  zero  except  win-m  -  0.  In  the  latter  case,  the  auto¬ 
correlation  r  (t  )  is  infinite,  hence  P  (t  )  can  be  represented  in  terms  of  the  Diroc 
Delta  function 

1'  (t  )  i'  •  Sfui  ■  &(i  )  (A  6) 
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where  k  is  a  constant  of  proportionality.  To  determine  the  value  of  k  and  to  show 
that  (A. 6)  is  the  correct  functional  form,  substitute  (A. 6)  Into  (A.3)orid  integrate; 

f 

S(u)  =  “  k  S(u)  •  /  6(t  )  •  cos  u  T  •  dt 

It  Jq 

or 

(A.  7) 

The  integral  in  tlie  denominator  of  (A  .7)  can  be  simplified  as  follows; 

6(t  )  •  cos  or  •  dT 

by  use  of  (A.2)above.  Thus  (A. 7)  becomes 

k  tr 

and  hence  the  autocorrelation  function,  1' (t  )  ,  for  white  noI';e  Is: 

r{T  )  IT  •  S(u)  •  6(t  )  (A. 3) 

Note  thnt 

SM  ^  S(f)  (A.  9) 

where  S(oj)  has  the  units  of  (psf)  /(rad/scc)  and  S(0  has  the  units  of  (psO'/cps. 
Substituting  (A. 9)  into  (A. 8)  gives  an  alternate  form  for  the  autocorrelation  function, 
namely, 


< 
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APPENDIX  B 


AUTOCORRELATION  FUNCTION  FOR  BANDWrPTH  LIMITED  WHITE  NOISE 

Given  o  random  process  which  is  characterized  by  bondwidth  limited  white  noise. 
The  power  spectral  density  (PSD)  for  such  a  piocess  is  denoted  by  S{u)}  and, as 
shown  in  Figure  19  S(u)  Is  defined  as 

S(u)  "  S  o,  /o  /u 

L  ^  ^  u 

-  S  .  -u  <u  <'-u, 
u  L 

=  0  otherwise 

S  =  constant  ^0 


From  the  WIencr-KInchIn  relations  (page  67  of  Reference  10)  the  autocorrelation 
function,  I  (t  ),  for  the  random  piocess  Is 


i  (t) 


du 


Ism 


u  T 

u 


sin 


(B.l) 


sin 


u 

u 


2 


T 


U 


cos 


T 
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or 


PM=  S  .  (w  -  u  )  » 

u  L 


■  .  “l 

sin 

1= 

•4 

_ 1 

u  +  u, 
u  L 

u 

2 

u 

L 

2 

A  graph  of  Tfr  )/S  .(u^  -  u^)  rs  shov/n  !n  Figure  20.  The  values  of  t  ^  and  t 
shown  in  this  graph  ore 


1  u  ^  o, 
u  L 


and 


2tt 


2  o  “  u, 

u  L 


The  following  properties  of  )  are  to  be  noted: 

1  .  Lim  I  (t  )  -  0 

T--»^ 


2. 


U),  /  w 
L  u 


1(0) 


S  •  (u,  -  u  ) 

L  u 


cos  Ut;  U  O,  ~  U 

U  L 


■  mean  square  value 

I'(t  )  j^S  •  (u^  -  U|^)j  • 

This  is  the  form  of  the  correlation  function  for  a  sinusoidal  signal;  and 

[S  •  (u  -  w,  )!  is  the  mean  squore  value  of  the  amplitude  of  the 
u  L 

signal . 

From  the  above  equation  for  t 


2 


as  u  — ►  (0,  (narrow  band) 

u  L 


T  2  — ♦O 


Gi  (u  -  u,  )-*-C3o  (broad  band) 
u  L 


(B.2) 


Also  for  T 


T  , — ►O  as  (u  ‘  u,  )  -*-00  (brood  borKl) 

I  u  L 


5.  Using  Equation  (B.1)  above,  it  Is  seen  that  when  =  0,  the 


autocorrelation  function  ('(t)  becomes 


1(1)'$  '  u 


sin  u  T 
u 


Figure  16  shows  two  graphs  of  i  (t  )  /  S  'u  ,  one  for  a  small  value  of  o  and  the 

u  u 

other  for  a  large  value  of  u  .  Negative  values  of  t  ore  shown  in  order  to 

u 

emphasl2e  the  symmetry  of  I  (t  ).  It  Is  seen  that  as  o  becomes  large  (large 
bandwidth),  the  correlation  function  decays  rapidly  with  time  delay  i  . 


For  Ihe  ('ase  of  Ihoorotical  while  noise. in  which  -  0  and 
-*■<»/  the  autocorrelation  function  approoches  the  Dirac  delto 
function,  6(t).  This  function  is  defined  as  follows:  (Reference  11) 


6(t  )  ■  Lir 


Si  n  af 


Thus  in  the  limit,  1'(t  )  becomes 


I  (t  )  ir  S  *6(7)  u,  -  0  ,  u  — ►GO 
_  w  L  u 


S  denotes  S  in  units  of  (  )^  /  rad  /  sec . 
o 
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Note  that 


Note: 


or 


S(f)  -  2n  SM 


S  -  2tt  S 
f  u 


2 

denotes  S  in  units  of  (  )  /  cps 


Thus  on  alternate  form  for  )  is 


r(T )  -  -j  '  ^  *‘’L  “ 


An  alternate  expression  for  the  Diroc  delta  function  can  be  obtained 
directly  from  the  above  Integral  expression  for  r(T  ): 

00 

r(T  )  -  s  .  f  cos  UT  ’  d(J 


“  I 


It  •  5  •  6(t  ) 

u 


or 


^7 


8(t  )  —  /  cos  U7  •  d(j 

■f) 


6(,  )  2] 


cos 


?TtfT  df 


or 


00 

L 


6(t  )  '  I  cos  2nfT  •  df  • 
'-00 


■rn. 
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Important  properties  of  8(t  )  ore 


Figure  2;  Onc-Dt'^rc'e-of-Frcedoni  Sprinp-Mass  SysttMii,  with  Fpriiig 

StiffiK'SS  k,  Danipinp  Constant  c.  Mass  m.  Applied  Force  F(t) 
and  Deflection  \V(t). 


F{»)  W(t) 


-♦ 

y 


Figure  3:  Uniform  B<'am  with  Point  Force  F(t)  Applied  at  Point  iTand 
Acceleration  Response  \V(t)  at  v. 


Ill 


F^(t)  W(0 


Figure  4:  Hif'id  I?ai'  Froo  to  TraiislaU*.  Init  CoiiKti  aincfl  in  Pitcli,  and 
Ilavinjf  a  I.ai'^;f'  Nnmlror  of  Appliod  Forca  s  F  |^(0 


F  ft) 

n 


Figure  5t  Uniform  Beam  Schematic  for  Response  at  ond'T'  on  A 

Excited  By  a  Point  Force  F  ft)  at  ~y  . 

n  '  n 
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Figurp  7.  Two  Single -degree-  of  Freedom  Systems  with  Impjise  Excitations 


Figure  8:  Graph  of  the  Complete  Crosscorrelation  Function,  Showing  the 
Slope  Discontinuity  at  r  -  0,  and  Showing  the  Possible  Different 
Forms  of  the  Crosscorrelation  for  Positive  and  Negative  Delay  Times 


jre  9:  Inputs  and  Outputs  For  the  Two  Single  Degree  of  Freedom  Systems 
Subject  to  Bandwidth  Limited  V/hite  Noise. 


0 


Frequency  (cps) 


Frequency  (cps) 


Frequency  (cps) 


Frequency  (cps) 


Figure  1 0; 


Illustrative  Exarrpie  of  the  Spectral  Density  S,  .(f)  and 

(r  -s) 

Tronsmission  loss  Coefficient  TL,  ,  of  tlie  Diffe  ence 

(i  -s) 

*  *  '  • 

Betweetr  Two  Responses  W  (t)  -  W  (t). 

r  s 
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in  rms  iounc!  Kressure 


(a)  Comparison  of  Slopes  of  Trend  Curves 


75  150  300  600  1700  2400 

Octove  Rands 


Roof  mean  square  acceleration  rotio 


Flgii'O  14  Lo  jncli  Acco  lorat  ions  Obtained  on  Structure  in  Forvvoid 

Half  of  Veil ic If  . 


•  □ 


foice  -fo  -  acceleration 
tranifer  function:  I 

nr 


Figure  17:  Dioginm  of  Seiies  Form  of  Force  -  to  -  Acceloiation 

~T 

Ttonsfei  Function  o  in  Terms  of  Acceptance, 
nr 

'rt  iinsniission  and  Rcrcivci'  riinctinns. 


122 


S(u) 


I  S 


-u 


“U, 


0  u 


u 


Figure  19:  Plot  of  Bondwidth  Limited  White  Noise  PSD 


Figure  2C  :  Graph  of  the  Autocorrelation  Function  for  Bondwidth 
Limited  White  Noise 
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